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SUMMARY 


This is the 4 th volume of a 5 volume report on the aerodynamics 
and aeroacoustics of advanced turboprops. It explains the use of 
the UAAP (Unified Aero-Acoustic Program) code for predicting 
overall propeller performance and for predicting steady and 
unsteady air loading, wakes, and noise of the blades. The 
aerodynamic and acoustic methods are based on linear pressure 
potential theory with corrections for non-linearity associated with 
axial mass flux induction and vortex lift on the blades. Theory is 
derived in Volume I (NASA CR-4329) and compared with experimental 
results in Volume III (NASA CR-185193) . Other volumes in the 
report are Volume II, which gives theory for shielding of propeller 
sources by a swept wing in compressible flow, and Volume V, which 
presents theory for acoustic shielding by a fuselage boundary 
layer. Theory in the latter 2 volumes is not included in the code 
described herein. 

Major sections of this volume include a general code overview, 
descriptions of the input and output, a general system description, 
and a listing of error codes. 
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SECTION 1 
INTRODUCTION 


This is Volume IV of a five volume report. The five volume report 
presents a unified theory for aerodynamics and noise of advanced 
turboprops. Aerodynamic topics include calculation of performance, 
blade load distribution, and non-uniform wake flow fields. Blade 
loading can be steady or unsteady due to fixed distortion, 
counter-rotating wakes, or blade vibration. The aerodynamic theory 
is based on the pressure potential method and is therefore 
basically linear. However, non-linear effects associated with 
finite axial induction and blade vortex flow are included via 
approximate methods. Acoustic topics include radiation of noise 
caused by blade thickness, steady loading (including vortex lift) , 
and unsteady loading. Shielding of the fuselage by its boundary 
layer and the wing are treated in separate analyses. 

The Unified Aero-Acoustic Program (UAAP) code discussed in this 
volume calculates the airloads on a single rotation Prop-Fan, or 
propeller, and couples these airloads with an acoustic radiation 
theory, to provide estimates of near-field or far-field noise 
levels. The steady airloads can also be used to calculate the 
non-uniform velocity components in the propeller wake. 

The airloads are calculated using a three dimensional 
compressible panel method which considers the effects of thin, 
cambered, multiple blades which may be highly swept. These 
airloads may be either steady or unsteady. The acoustic model 
uses the blade thickness distribution and the steady or unsteady 
aerodynamic loads to calculate the acoustic radiation. 

This report constitutes the users manual for the UAAP code. The 
report is divided into five sections: 

General Code Description 
Input Description 
Output Description 
System Description 
Error Codes 

The user must have access to IMSL10 libraries (MATH and SFUN) for 
numerous calls made for bessel functions and matrix inversion. 

For plotted output users must modify the "dummy" calls to 
plotting routines included in the code to system-specific calls 
appropriate to the user's installation. 

Instructions for use of the separate code for prediction of 
boundary layer refraction are included in Volume V which 
describes that code. 
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SECTION 2 

GENERAL CODE DESCRIPTION 

The UAAP code is divided into six major functional modules. Each 
functional module performs a unique task and consists of one or 
more subroutines. Generally, the routines in a functional module 
are unique to that module although, interpolation routines, 
integration routines and others may be required in another 
functional area. The functions performed by these six modules 
are: 

1) Blade Geometry Generation - Inputs the propeller/blade 
shapes and transcribes them into the format required by 
the other portions of the code. 

2) Propeller Flow Field - Inputs the propeller axial 
velocity flowfield (supplied by the user) and transcribes 
that flowfield into the format required by other portions 
of the code; interpolates the propeller flowfield at the 
control points required by the aerodynamic portion of the 
code. 

3) Compressible Panel Method Aerodynamics - Calculates the 
steady or unsteady aerodynamic loads on the blade. 

4) Vortex Flow Aerodynamics - Calculates the change in the 
airloads due to leading edge and tip edge vortex effects 
produced by the steady compressible panel method. 

5) Acoustic Analysis - Calculates the near field or far 
field propeller noise. 

6) Wake Calculations - Calculates the wake caused by viscous 
and potential effects. 


Options have been incorporated in the code to allow restart of 
the code at two points: 

1) Restart of Compressible Panel Method Aerodynamics 

2) Restart of Acoustic Analysis 


The restart capabilities are detailed in the Input Description 
section of this manual. 

This code involves contributions from several disciplines in 
propeller design, including performance, acoustics, and 
structures. Each of these disciplines has its own coordinate 
system and these are supported whenever possible. However, this 
does bring about some conflict and the user is cautioned to 
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reference this users manual for proper interpretation of the 
coordinate systems. 

Figure 1 presents an overview of the UAAP code. As shown, two 
parallel paths exist within the code: a steady airload 

computation, and an unsteady airload computation. These 
computations can only be performed with separate executions of 
the code. 

The steady airload calculation path shows the computation of the 
loading kernel matrix and thickness vector, and the storage of 
those quantities in an external file. The stored loading kernel 
and thickness vector are used to compute the blade loading 
including the vortex loading. 2-D airfoil tables are used to 
obtain the viscous losses. The blade loading is then used to 
obtain blade wake profiles, overall aerodynamic performance, and 
noise. 

The unsteady airload calculation follows a path similar to that 
for calculation of the steady airloads. As shown in figure 1, 
the unsteady loading does not include vortex loading or viscous 
drag, and produces only an unsteady acoustic calculation. Complex 
values of the inlet downwash vector are required as input for the 
unsteady calculation. 

Steady and unsteady airload and acoustic calculations must be run 
separately, thus, output from the unsteady calculation doesn't 
contain the steady loading effect. Also, unsteady loads must be 
run one harmonic at a time. 
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SECTION 3 
INPUT DESCRIPTION 

The primary input mode is expected to be a predefined input 
dataset. Upon program execution the primary input file, file #5, 
is immediately read and echoed to both the primary output file, 
file #6, (to obtain a complete copy of the input on the output 
file) and to a scratch file, file #11. All subsequent primary 
input is obtained by reading file #11. 

No other input files are required by the code except as the user 
may define for restart capabilities. Several scratch files are 
generated by the program? they are defined in the section "UAAP 
System Description." 

There are many input values which act as tolerances on 
iterations, or convergence criteria. These have been defaulted 
within the code to "recommended" values. Although they 
explicitly appear in the input description, they need not be 
input. Input data which need not be input appears at the end of 
each section, and the input location number appears within 
parentheses. 

The next sections present a description of all the input 
necessary to run the code. 
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Commands 


The code is command driven. Commands are used to identify the 
input sections, the input, and to execute various functions of 
the code. Sub-commands are used to execute options within 
functions. The general input format of the command is: 

(starting in column 1) 

COMMANDX ( SUBCOMMD) where "COMMANDX" is an 8 character word, 

including trailing blanks, and "SUBCOMMD" 
is also an 8 character word, including 
trailing blanks. 

Load 

Some commands require that numerical input must be read in next. 
For this purpose a location-specified free field input routine is 
used which reads between columns 1 and 72. Locations for input 
are indicated by "L" followed by a number indicating the desired 
start of a location. 

Some of the location specified input controls integration mesh 
sizes, Fourier series convergence, and other program tolerances. 
These locations have been identified by parenthesis () around the 
location number and can usually be ignored by the user. However, 
if these values are input, the parenthesis () should be omitted 
from the location field. 

Input for this "LOAD" format is illustrated below: 

(starting in column 1) 

C THIS IS AN OPTIONAL COMMENT RECORD (or Records) 

L 10 1.1 2 3.5 

C THIS IS ANOTHER OPTIONAL COMMENT RECORD (or Records) 

12.7 13.1 

L 91 6.5 

END 

This will cause the locations shown below to have the following 
values : 


D CATION 

VALUE 

10 

1.1 

11 

2.0 

12 

3.5 

13 

12.7 

14 

13.1 

91 

6.5 


Notes: "Scientific notation" is not allowed, e.g. 2.1 E+03 is 

not allowed. Implied repetitions are not allowed, e.g. 3*2.1 is 
not allowed. The "END" record is required to terminate each 
entry into the LOAD routine. 
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Command Summary 


A list of accepted commands are shown below. Since this is a 
command driven code, command order is important, and therefore 
the commands are listed below in the required logical sequence 
for execution. 


HEADER 

RUNPARMS 

AIRPARMS 

BLADEGEO 

LS TP ARMS 
NOIZPARM 

VELGRADS 

VORTPARM 
WAKE P ARM 
AEROEXEC 
WAKEEXEC 
NOIZEXEC 
ENDCASE 

ENDJOB 


Input page header cards 

Input flight parameters 

Input options for 2-D drag look-up. 

Input propeller/blade geometry (2-D,RXY or 
XYZ coordinates) . 

Input options to panel aero code. 

Input options to the noise portion of the 
code. 

Input axial velocity field on a defined 
grid. 

Input options to vortex flow calculation. 

Input options to wake calculation. 

Execute the aero code. 

Execute the wake calculation. 

Execute the noise portion of the code. 

End of input and calculations for the 
current case. 

End of job, terminate the program. 
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Specific Input Description 


Input requirements/ options are provided below for each command, 
in the order listed above. 


HEADER : 

This command enters page heading records. Records following the 
HEADER command are sequentially read until an "END" command is 
found. Up to 10 page heading records may be entered. 


RUNPARMS : 

This command enters run parameters using "LOAD" format. The 


location, 

default, 

and input 

descriptions are: 

LOC. 

Default 

Variable 

Description 

1 

0 

RUNBUG 

Debug Option, 1 turns on. 

2 

59.0 

QDEGF 

Ambient temperature, °F. 

3 

1.0 

QRHOR 

Ambient density/Sea Level Std. 
density. 

4 

0 

QADVF 

Advance ratio, J = V^nD. 

5 

0 

QMX 

Free stream Mach number. 


LSTPARMS : 

This command enters input options to the panel aero portion of 
the code using "LOAD" format. 

The location, default value and input description are: 

( * See Figure 2 for pictorial description) 


LOC. Default Variable Description 


1 

* 2 


0 

10 


ZSTBUG 

QNCP 


8 


QNSM 


20 3 QPART1 


Debug Option, 1 turns on. 

Number of chordwise panels 
(maximum of 20) 

Number of spanwise modes 
(maximum of 20) 

0 generate and W T 

1 generate K~ and read W T 

2 read k” and gen. W T 

3 read in k” and U T 

The code requires an inverse 
kernel matrix (K~ ) and a 
thickness vector (W T ) to obtain 
the aerodynamic loading on the 
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21 


22 


28 


29 


* 101 


* 141 


181 


351 


0 QQ 

0 QK 

2 QNBOPT 


10 QITNON 


.2, .35, QZAR 

.45, .55, 

.65, .75, 

.85, .95 

. 5 , . 5 QCONTP 


. 4 , . 4 QCHW 


, 0 QWMU 


blade^. 

The K matrix and w^ vector 
require a significant amount 
of CPU time, and an option to 
utilize a previously generated 
K matrix and W T vector has 
been incorporated into the code. 

See SECTION V. 

Order of unsteady loading 
harmonic, 

= 0 for steady loading. 

Number of circumferential modes 
for unsteady loading, = 0 for 
steady loading 

= 1 for supersonic leading edge 
element, use when supersonic flow 
expected at leading edge. 

= 2 for subsonic leading edge 
element 

Max number of non-linear 
iterations. 

=0 for linear calculation. 

Spanwise locations of control 
point radius. 

The number of these must agree 
with QNSM, fraction of RTIP. 

Chordwise location of control 
points within each panel 
normalized to the panel width. 
There must be QNSM of these. 

Width for chordwise averaging of 
downwash, normalized to panel 
width. There must be QNSM of 
these input. 

Downwash vector for unsteady 
loading calculation. These can 
be obtained from Appendix I, 
equation (8) . 

The values of QWMU must be input 
sequentially starting in 
location 351, as real part, 
imaginary part for each control 
point across the 1st spanwise 
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station, and then proceeding 
outward along the blade span. 
Thus, there must be (QNCP*QNSM) 
pairs of values input for QWMU. 

The following input locations may be ignored. 


(3) 

.02 

QDR 

Radial extent of singular 
integral at control point radius 
for wake calculation radius. 

(5) 

1024 

QN 

Number of points in FFT for 
terms in kernel integration, 
max. of 2048, must be a power of 
2. 

(6) 

.004 

QDELTA 

Axial step size for FFT. 

(7) 

4 

QMODOP 

Spanwise mode shape option 
(use default value) . 

(8) 

.0001 

QTOLF 

Tolerance for W (omega) 
integration. 

(9) 

.001 

QTOLT 

Tolerance for other 
integrations . 

(10) 

.005 

QTOLS 

Tolerance for summations 

(11) 

20 

QMM1 

Loop limit for harmonic 
summation in wake kernel. 

(12) 

30 

QMM2 

Loop limit for harmonic 
summation in bound kernel. 

(13) 

10 

QMM3 

Loop limit for harmonic 
summation in thickness vector. 

(14) 

8 

QITABK 

File number allocated for K-l 
storage . 

(15) 

9 

QITABT 

File number allocated for W T 
(thickness vector) storage. 

(16) 


QINT4 

Not currently used. 

(17) 

0 

QPRINT 

Additional debug print - not 
recommended . 

(18) 

0 

QPRIN1 

Additional debug print - not 
recommended . 
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(19) 

0 

QPRIN2 

Additional debug print - not 
recommended . 

(23) 

.025 

QKDOWN 

Radial step size for 
non-singular bound kernel 

(24) 

.010 

QKSTART 

Width of singular region for 
bound kernel integration 

(25) 

5 

QMM4 

Loop limit for summation in 
sound power calculation. 

(26) 

.9 

QMBLEN 

Trailing edge effective Mach 
number for blending to 
supersonic trailing edge 
elements. 

(27) 

1024 

QNO 

Number of points in FFT 
integration for n = 0 term in 
kernel, max =2048. 

(221) 

.002 

QINMES 

Step size for radial integration 
in wake kernel. 
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BLADEGEO: 

This command enters the propeller/blade geometry using the "LOAD" 
format. For any propeller and flight condition, the shaft power 
required, and the thrust produced by the propeller are a function 
of the blade angle. In using this code it is recommended that the 
blade angle be adjusted so that the power or thrust calculated by 
this program matches a desired value, (see Volume III, page 6.2 
for further explanation) . 

Three subcommands exist to allow input of the geometry in 
different forms: 

BLADEGEO ( 2 -DCOORD) , BLADEGEO (RXYCOORD) , and BLADEGEO (XYZ COORD ) . 
These are described below. 

BLADEGEO f 2 -DCOORD) : 

This form of the blade geometry is expected to be the most widely 
used. Input takes the form of the spanwise variation of 
thickness ratio, chord/diameter ratio, twist, airfoil section 
designation and stacking axis coordinates. This command 
(sub-command) will calculate the blade surface coordinates and 
interpolate these into a form required by the other parts of the 
program. To assure that there are no errors in the blade output 
description only a limited amount of extrapolation of the input 
blade coordinates is allowed. Thus, it is best to provide input 
stations, X, and streamline angles, SLA, which will define the 
root and tip sections of the blade such that the output stations, 
ZBLDST, can be interpolated and not extrapolated. The code 
requires exactly ten (10) inputs defining X. The input described 
below is illustrated in Figures 3 and 4. 

Note that there are empty locations in the input. These are not 
used by the program. 


LOC. 

Default 

Variable 

Description 

1 

0 

BLDBUG 

Debug option, 0 is off, 1 is 

31 

0 

BLADE 

Number of blades 

32 

0 

D 

Propeller diameter, ft. 

33 

0 

SCO 

Propeller hub/tip ratio 

41-50 

0 

X 

Spanwise input stations, 10 
required, fractions of RTIP. 

51-60 

0 

HOB 

Spanwise airfoil maximum 
thickness/chord ratio. 

61-70 

0 

BOD 

Spanwise chord/diameter 
variation. 
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71-80 


0 


CLD 

81-90 0 DTHET 

347 0 THTDES 

348 0 THTCUT 


381-390 0 SLA 


711-720 0 XSWP 


721-730 0 YSWP 


731-740 0 ZSWP 


Spanwise variation of design 
lift coefficients. 

Spanwise twist variation. 

Twist should be input such that 
the 75% radius twist is 0, 
degrees . 

The 75% radius value of blade 
angle at which the blade is 
to be "designed", degrees. 

The 75% radius value of blade 
angle at which this calculation 
is to be run, degrees. 

Note: Propeller blades are 

assumed to be "designed" at one 
value of blade angle. The input 
values of thickness/ chord , 
chord/diameter, camber, twist 
and stacking line are assumed to 
be defined at the value of 
blade angle in Location 347. 

Location 348 defines the blade 
angle at which this calculation 
is to be run. 

Spanwise variation of streamline 
angle, degrees. Blade airfoil 
sections are assumed to be on 
cones which approximate the 
streamlines thru the propeller. 
This input is the 1/2 cone 
angle, positive as shown in 
Figure 3 . 

Spanwise variation of X 
coordinate of mid-chord stacking 
line, fraction of Rtip, see 
figure 4. 

Spanwise variation of Y 
coordinate of mid-chord stacking 
line, fraction of Rtip, see 
figure 4. 

Spanwise variaton of Z 
coordinate of mid-chord stacking 
line, fraction of Rtip, see 
figure 4. 
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900 

21 


ZNPCOV 

Number of output stations in 
the chordwise direction, max 49 


901-949 

0 . 

, 5. 

, PCTCHD 

Chordwise location of output 



10, 

15, 

20, 

stations. Both upper and lower 



25, 

30, 

35, 

airfoil surfaces are output at 



40, 

45, 

50, 

the same chordwise locations, 



55, 

60, 

65, 

% chord. 



70, 

75, 

80, 


— 


85, 

90, 

95, 




100 




— 

950 

21 


ZNIS 

Number of output stations 
in the spanwise direction, max 
of 49. 

— 

951-999 

.10 

, .30 

, ZBLDST 

Spanwise location stations at 


.40, .50, which the blade surface will be 

.55, .60, defined, radius ratios, r/RTIP. 

.65, .70, 

.75, .80, 

. 825 , . 85 , 

.875, .90, 

.925, .95, 

.96, .97, 

.98, .99 

1.0 

The following input locations may be ignored. 

(346) 2 SWPOPT Input sweep option. The default 

value is required. 

(349) 1 ZKCUT Type of airfoil section defined 

0 planar, 

1 conical, use default, 

2 cylindrical 

(741-751) 1 BAFL Integer characterizing the 2-D 

airfoil at each spanwise 
station. 

Only NACA 16 series airfoils can 
be generated with this deck. 
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BIADEGEOfRXYCOORm : 

An optional method for input of the blade geometry has been > 
provided. This form allows for up to 50 input stations. 
Additionally, it requires a table of the displacement of the mean 
camber line from the chord line. 

This option requires that the blade sections being described be 
on cylinders whose axis is the centerline of rotation. The input 
is illustrated in Figures 5, 6 and 7. 


LOC. 

1 


Default 

0 


4 

5 

6 


0 

0 

0 


51 -100 0 
101-150 0 


Variable 

BLSBUG 

STANO 

PCTNO 

DIAMET 

SPINNER 

BLADES 

CTSTA 

PCTCD 


At each spanwise station 


151-200 0 


201-250 0 


251-300 0 


THROB 


CHDOD 


CAMBR 


301-350 0 
351-400 0 

401-450 0 


TWIST 

XMC 

YMC 


Description 

Debug Option, 0 is off, 1 is on. 

Number of spanwise input 
stations, max of 50 

Number of chordwise input 
stations, max of 50 

Propeller diameter, ft. 

Propeller hub/tip ratio 

Number of propeller blades 

Spanwise input stations, r/R TIp 

Chordwise input stations, %chord 


Maximum blade thickness, 
fraction of chord. 

Blade chord, fraction of 
diameter. 

Equivalent NACA Series 16 
camber. Obtained by taking the 
non-dimensional maximum height 
of the blade mean camber 
line / .05515 , fraction of 
chord . 

Blade section chord angle, deg. 

X-coordinate of mid-chord 
stacking line, fraction of Rtip. 

Y-coordinate of mid-chord 
stacking line, fraction of Rtip. 
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451-500 

0 

ZMC 

501-550 

0 

XSLA 


1001 1 
1002 1 


Z-coordinate of mid-chord 
stacking line, fraction of Rtip. 

Streamline angle, deg. 

This is the flow angle relative 
to the centerline of rotation, 
and is used in the sweep angle 
calculation. 

Required value, to set up table 
for interpolation. 

Required value, to set up table 
for interpolation. 


Blade mean camber line displacement table. 


1003 0 


1004 0 


1005 0 


1005 0 

+ 

PCTND 


1005 0 

+ 

PCTND+STAND 


PCTND Number of chordwise points in 

mean camber line displacement 
table. 

STAND Number of spanwise points in 

mean camber line displacement 
table. 


CAMLN Ascending array of chord 

fractions, PCTND values. 

Followed immediately by : 

CAMLN Ascending array of spanwise 

radius, fraction of Rtip. 
Starting after last input 
location, STAND values. 

Followed immediately by : 


CAMLN (I,K) 
1=1 , STAND 
K=1 , PCTND 


Array (PCTND*STAND) of 
mean camber line displacements 
from chord line at each radial 
station, at 1st chordwise 
location, fraction of chord. 
This is followed by a 
similar array for the 2nd 
chordwise location, and 
continues thru the PCTND 
chordwise location. 
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BLADEGEOfXYZ COORDS : 

This form of blade geometry inputs the blade description in XYZ 
coordinates of each blade surface. Figure 8 illustrates the 
coordinate system used. Note that this XYZ coordinate system is 
different than the 2-D or the RXY system. 


LOC. 

Default 

Variable 

Description 

1 

0 

BLDBUG 

Print option, 0 is minimal 
1 is additional. 

2 

0 

DIAME 

Propeller Diameter, ft. 

3 

0 

SPINN 

Hub to tip ratio. 

4 

0 

BLADE 

Number of blades 

5 

0 

STINN 

Number of radial stations in 
input. 

6 

0 

CDINN 

Number of chordwise stations in 
input . 

7 

0 

STOUTN 

Number of output radial stations 
needed to define the blade. 

8 

0 

CDOUTN 

Number of output chordwise 
stations needed to define the 
blade. 

51 

0 

STOUTV 

Values of radial stations, 
fraction of radius. 

101 

0 

CDOUTV 

Values of chordwise stations 
fraction of chord. 


An "END" record is required to terminate the "LOAD" input. The 
"END" record is immediately followed by the xyz coordinates of 
the blade as follows: 

For each of the STINN radial input stations the following 
input records are required: 

a) A label record 

b) For each of the CDINN chordwise input stations the 
X, Y, & Z coordinates of the face (pressure) side of 
the blade are input. With one set of X, Y, Z, 
coordinates in free field format per record. The 
units are inches. 
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c) Another label record 

d) For each of the CDINN chordwise input 
stations the X, Y, & Z coordinates of the camber 

(suction) side of the blade surface are input. 

With one set of X, Y, Z coordinates in free field 
format per record. The units are inches. 
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VELGRADS : 

This command will input or initialize the axial velocity ratio at 
the propeller. These ratio may be due to the spinner-hub and or 
nacelle effects on the freestream flow, but the result must 
appear to the code as an axisymmetric flowfield. 


LOC. Default Variable 

1 0 V1BUG 

2 0 V10PT 

For option V10PT=1. only: 

3 0 VINRD 

4 0 VINAS 

5 0 VIVRF 

6 0 VITIP 

7 0 VICLL 


Description 
Debug option 

Velocity ratio input option, 

0 : initialize velocity field 
to 1.0, fraction of freestream 
velocity. No further input is 
required. 

1 : input the velocity ratio. 

The axisymmetric velocity ratio 
is input on a grid of radial and 
axial points which encompass the 
blade outer boundaries. The 
flowfield values are 
interpolated at specific "nodal” 
points required by the code. 

Number of radial stations in 
grid. 

Number of axial stations in 
grid. 

Reference velocity, by which 
local velocities are normalized, 
same units as VIRD, VI W, and 
VIAX. 

Reference radius, by which local 
radii are normalized, same 
units as VIRD and VIAX. 

Centerline location. The axial 
location of the blade pitch axis 
is assumed to be 0. This value 
will be subtracted from the 
input values of VIAX after 
input so that the axial location 
of the blade centerline will be 
at 0. Same units as VIAX. 
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The following inputs are 
repeated for each of the axial 


stations in 

the grid, (K=l 

, VINAS) . 

25 0 

VIRD (I , K) 
1=1, VINRD 

Radial locations of 1st set of 
up-stream grid points, 
same units as VITIP. 

25+ 

VIW ( I , K) 
1=1, VINRD 

Velocities at VIAX (I) , VIRD (I) , 
same units as VIVRF. 

25+ 

VIAX (I , K) 
1=1, VINRD 

Axial locations of VIW (I), 
same units as VITIP. 


WAKE P ARM : 

This command provides input parameters to the propeller wake 
calculation procedure. This input is used in both the 
performance calculations, and in the "Wake" calculation option. 
The variables used for propeller efficiency and wake calculation 
are listed below. 


LOC. Default Variable 

1 0 WAKBUG 

6 1 YIPLOT 

8 1 YNPSKN 

20 0 YIPRT 

23 0 YIVWK 

25 11 YNORPE 


Description 

Debug option, 1 print 
WAKEPARM input data. 

0 do not print data. 

Wake output, 1 is on, 0 
is off. 

Skin friction drag . 

1 : included, 

0 : not included. 

Output option: 

0 : no output 

1 : Vr/V & Vx/V 

2 : Fluctuating lift 

The viscous wake shape is 
described by 1 of 3 
options 

0 : gauss pulse 

1 : cosine squared 

2 : cosine 

No. of output radii. 
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751 


Output radii, r/Rtip. 


.24, .35, YROUTP 
.45, .55, 

. 65 , . 75 , 

.80, .85, 

.90, .95, 

.98 


The following input locations may be ignored. 


(2) 

.05 

YK 

Spanwise integration mesh 
size, fraction of radius. 

( 9 ) 

.0001 

YTOL 

Fourier series sum 
convergence tolerance 

(10) 

99 

YMM1 

Max. no. of Fourier 
coeff. to calculate, max. 
of 99. 

(12) 

.001 

YTOL1 

Tolerance for high 
frequency form of Fourier 
coeff. 

(13) 

0 

YLASTM 

Last Fourier coeff. for 
detailed output. 

(16) 

0.7 

YZO 

Origin of special routine 
for interpolation of 
circulation curve. 

(17) 

0.5 

YZNORM 

Circulation curve 
normalizing factor 

(18) 

1.0 

YAOPT 

0 to not iterate on 
induced angle 

(21) 

0 

YMM2 

Turns on diagnostic print 


for Fourier coeff. up to 
YMM2 when YIPRT = 2. 
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The following inputs are required in addition to the above if 
wakes are to be calculated. 


LOC. 


4 


5 


50 


51 


52 


53 


351 


701 


711 


Default Variable 


0 YNOR 

0 YNPX 

0 YXNCS 

0 YXMMU 

0 YXMMV 

0 YXMMW 

ROUT 

0 YARRY (I , K) 

1=1, YXNCS 

0 YARRY (I, K) 

1=1, YXNCS 
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Description 

Output radii are divided 
into those within the tip 
radius, YNOR, and those 
beyond the tip radius, 
YNPX. 


Number of output radii 
which are less than or 
equal 1.0 (within tip 
radius) . 

Number of output radii 
which are greater than 
1.0 (outside of tip) . 

No. of output axial 
locations, max of 10 at 
each radius. 

No. of Fourier coeff. 
used in calculation of 
chordwise wake component. 

No. of Fourier coeff. 
used in calculation of 
radial wake component. 

No. of Fourier coeff. 
used in calculation of 
downwash wake component. 

Output radii, fraction of 
radius. There should be 
YNOR+YNPX of these. 

Output axial locations, 
fraction of radius, for 
1st output radius. There 
should be YXNCS of these. 

Output axial locations, 
fraction of radius, for 



2nd output radius. There 
should be YXNCS of these. 


721 0 


YARRY (I , K) Output axial locations, 

1=1, YXNCS fraction of radius, for 

3rd output radius. There 
should be YXNCS of these. 


where 

The following input 
(54) 0.01 


K = 1, YNOR + YNPX 

location may be ignored. 

YFPHW Width of averaging 

function for wake 
calculations = X/R TI 


AIRPARMS: 

This command will input options which are used to obtain the 2-D 
airfoil drag from built-in tables of lift and drag coefficients. 


LOC. 

Default 

Variable 

Description 

1 

0 

AIRBUG 

Debug Option, 1 to print 

The following input 

locations may be 

ignored . 

(2) 

24 

AIRNUM 

2-D Airfoil data pack 
number wanted for this 
run. 

Currently, only NACA 16 
is available in this 
code. 

(3) 

5 

AIRTYP 

0 : Isolated 2D airfoil 

1 : Isolated 2D airfoil 
w/cascade corr. 

4 : Isolated 2D airfoil 
w/ sweep corr. 

5 : Isolated 2D airfoil 
w/ sweep & cascade. 

(13) 

0 

AIRCOR 

0 : no drag correction 
input 

1 : drag correction 
input 

(15) 

0 

CDMLT 

2D drag multiplier, used 
if AIRCOR = 1. 

(17) 

0 

DCDT 

24 

2D drag increment, used 
if AIRCOR = 1. 



CDMLT and DCDT are used to alter the drag of the 
stored airfoil data. They may be used to simulate 
another airfoil, to account for a rough airfoil or 
provide a better match with test efficiency. 

C Q = C D * CDMLT + DCDT 
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NOIZPARM : 

This command will input the options which are used to control the 
acoustic calculations. 


NOTES: B is number of blades, M is noise harmonic order 

B * ZMMAX must be less than 1001 

MX is flight Mach number (may not be 0 for near field 
calculation) 

LOC. Default Variable Description 

1 0 ZNZDBG Debug option : 

0 no print; 

1 print input parameters 

2 print file 50 input 


2 0 ZNFIND 

3 1 ZMMAX 


0 for far field theory 

1 for near field theory 

Max noise harmonic order 
to calculate. 


4 ALT 

9 ZNX 


10 1 ZX0RX1 


Distance from prop axis 
to observer, ft. 

Number of axial 
directivities for 
noise calculations (20 
max) 

Axial directivities are 
1 : visual or 
0 : retarded . 


11 XI Visual positions of 

observer along axis 
( + ahead of prop), ft. 
(calculated if X is 
input) . 
or 

31 X Retarded positions of 

observer, ft. 
(calculated if XI is 
input) . 
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59 

1 

ZIBLW 

1 to add boundary layer 
and wake displacement 
to geometry of blade 
trailing edge. 

0 for comparison with 
other predictions. 

70 

0 

PHI 

Azimuthal observer 
directivity angle. 0 is 
when the blade lies 
in the plane of the 
observer and prop axis 

PHI - 

Azimuthal observer directivity angle for unsteady 
loading, degrees, PHI increases in the direction of 
rotation, zero corresponds to the phase reference 
position for calculation of unsteady airloads. Thus, 
the observer position is specified relative to the 
unsteady flow field reference. 

The following input 

locations 

may be ignored. 

(51) 

10 

ZJMINI 

With ZNJ1, controls 
chordwise integration 
mesh 

(52) 

8 

ZNJ1 

# chordwise points = 
ZJMINI + M + B * ZNJ1, 
1001 max pts 

(53) 

10 

ZJMIN2 

Used with ZNJ2, controls 
spanwise integration mesh 

(54) 

10 

ZNJ2 

# spanwise points is JMAX 
(must be less than 2001) 
which is calculated : 
if ZNJ2 > 0 

JMAX = ZJMIN2+M*B*NJ2 ; 
if ZNJ2 < 0, 

JMAX = ZJMIN2 - ZNJ2 
if ZNJ2 = 0, 

JMAX = ZJMIN2 + ZNJ2 
points/phase cycle. 

(55) 

6 

ZIW 

Output file # 

(56) 

1 

Z START 

Start harmonic order 

(57) 

1 

ZMINC 

Increment in harmonic 


order 
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(58) 


2 


ZLHSOR Indicates what type of 

unsteady flow 

1 = blade wakes, 

2 = flowfield, 


(60) 

0 

ZPRUNS 

If not 0 print diagnostic 
unsteady loading noise 
table 

(61) 

0 

ZITRAP 

If not 0 max diagnostic 
print 

(62) 

0 

ZNFOUT 

If not 0 print harmonic 
table to file # ZNFOUT, 
allocate appropriate 

(64) 

0 

ZKXTND 

See note below. 

(65) 

10 

ZKMIN 

See note below. 

(66) 

10 

Z OMEGA 

See note below. 

(67) 

0 

ZIPDET 

If not 0, print tables of 
amplitude and phase vs. 
radius. ZIPDET is the 
highest order to be 
printed. 

(68) 

0 

ZLHRMP 

If not 0, print load 
harmonic summation terms 
(far-field option only) . 

(71) 

1 

ZNAIR 

Code for airfoil 
thickness 

distribution (specify for 
all radial stations) : 

1 = series 16 


3 = series 64 
5 = series 65 
7 = biconvex parabolic 
(analytic) 

8=4 digit series 

ZKXTND, ZKHIN, and Z0ME6A s 

In the near field option, the program 
evaluates a Fourier transform numerically using 
rectangular integration. The range for this 
frequency integration is 1/(1+MX) to 1/(1-MX), where 
MX is the flight Mach number. The number of steps in 
the integration range determines the trade-off 
between precision and running time. The integrand is 
typically a fluctuating quantity whose rate of 
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oscillation dl/dw, is computed by the program as a 
function of several factors such as harmonic order, 
Mach number, and observer position. To achieve 
uniform precision over a range of conditions, the 
program determines the number of points, kk, in the 
frequency integration range from 2 input numbers: 
the first is the minimum points in the range (ZKMIN) 
and the second is the number of mesh points per 
oscillation of the integrand (ZOMEGA) . The number of 
integration points is then kk = ZKMIN + ZKOMEGA * 
dl/dw. Default values are ZKMIN = 10 and ZOMEGA = 

10, which numerical tests indicate gave a reasonable 
compromise between precision and running time. If 
the user wishes to experiment, he can override 
either of the 2 default values. Reducing the number 
of points will reduce running time and storage 
requirements, but also reduce the precision of the 
calculation. Determining the satisfactory level of 
tradeoff between precision and running time for 
other than the defaults is up to the user. 

In some cases a significant contribution to the 
noise can be caused by frequencies outside the range 
of frequency integration noted above. To account for 
this, ZKXTND points are included outside both the 
upper and lower bounds of the integration. If ZKXTND 
equals zero then it is automatically computed by the 
program, otherwise the input value of ZKXTND is 
used. However, kk + 2*ZKXTND must be less than 401. 
If this is not true then the code will reduce ZMMAX 
by one until this condition is met. 
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VORTPARM : 

This command inputs the options which are available to control 
the Vortex Flow Aerodynamic calculations. 


LOC. 

1 


2 


3 


4 


5 

6 


Default 

0 

1 

1 

1 


.97 

0 


Variable 

VTXDBG 

DO LEAD 
DOS IDE 

DOAUGL 

ZAUGFF 

TIPLOR 


Description 

Debug option 

if = 0 no input print, 

if - 1 input printed 

If = 1 then calculate 
additional lift due to 
leading edge vortex. 

If * 1 then calculate 
additional lift or 
radial force due to tip 
edge flow. 

If = 1 then calculate 
additional lift due to 
leading edge vortex 
shed over aft portion of 
blade at the tip 
(augumented lift) . 

Radius at which augmented 
lift acts, r/Rtip. 

Indicates whether tip 
edge flow results in 
extra lift (0) or radial 
force (1). 


TIPLOR is set by the user to determine the type of 
edge flow. 

0 for separated tip flow, where the tip 
vortex gives extra flow at the tip; 

1 for attached tip flow which produces a radial tip 
edge force. 
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AEROEXEC f PRNTCASE ) : 

This command (Subcommand) will process all of the input and 
execute the code to the point where the Compressible Panel Method 
prints the data it has received. The aerodynamic calculations 
will not be performed. This is useful in verifying the input 
before doing the actual aero calculations. 

AEROEXEC (EXECCASE1 ; 

This command (Subcommand) will processs all of the input and 
execute the aerodynamic calculations. 

NOIZEXEC f EXECCASE^ : 

This command will execute the noise calculation section of the 
code. This should be placed after the AEROEXEC record. 

WAKEEXEC ( EXECCASE ) : 

This command will execute the potential and viscous wake 
calculations. This should be placed after the AEROEXEC or 
NOIZEXEC record. 

ENDCASE : 

This command signifies the end of all commands for the current 
case, and is required. 

ENDJOB : 

This command will terminate execution, and is required. 
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SECTION 4 

OUTPUT DESCRIPTION 

This section presents a description of the output from the UAAP 
program and includes output print from AEROEXEC (the panel code 
performance calculation) , WAKEEXEC (the wake velocity 
calculation) , and NOIZEXEC (the noise calculation) . 

The debug option described under the various input section 
commands controls the scope of output, and are either 0. (the 
default), 1., or 2., yielding additional output. Most pages 
have a page title followed by the label records entered in the 
input HEADER section. The page title consists of : 1) a brief 
description of the page contents. This is followed by, after the 
first colon (;), the subroutine name which is doing the printing 
of this page, 2) the input command option currently in effect, 

3) the time and date at the start of the run, and 4) the 
program/version identification. 

Output Description ( DEBUG = l) 

OUTPUT ; 

The output starts with an echo of the input data set, which is 
discussed in the input description, and shown in Figures 9-13. 

RUNPARMS : 

The debug option of 1. will print the page shown in Figure 14, 
which shows the input parameters, defined in the input section, 
along with the location number and the input or defaulted value. 

NOIZPARM : 

The debug option of 1. will print the page shown in Figure 15 
which shows the input parameters, defined in the input section. 

VORTPARM : 

The debug option of 1. will print the page shown in Figure 16 
which shows the input parameters, defined in the input section. 

LSTPARMS : 

Selecting the debug option, LSTBUG = 1., will print additional 
input definition as shown in Figure 17, which lists the available 
input options, their location number and the selected or 
defaulted values. Refer to the input description for more input 
detail. 

AIRPARMS ; 

Figure 18 presents the additional output for debug = 1. Again, 
the parameters are defined in the input section. A few lines of 
airfoil description are given, defining the airfoil selected, 
used to determine profile drag. 
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WAKEPARM : 

If a debug option of 1. is chosen, an additional input table. 
Figure 19 is obtained. Refer to the input section for parameter 
definition. 

VELGRADS: 

Figure 20 is a sample of output using debug = 1 , which labels the 
input for this section. 


BLADE GEOMETRY OUTPUT 


BLADEGEO : 

DEBUG =1. in the BLADEGEO section of input yields an additional 
page of input definition, given in Figure 21. Definitions of the 
parameters follow: 


Parameter 

DEBUG 

BLADN 

DIAMETER 

SCO 

SWEEP TYPE 
DESIGN ANGLE 
RUNNING ANGLE 
TYPE CUT 
X 

T/B 

B/D 

CAMBER 
DELTA THETA 

XSWP 

YSWP 

ZSWP 


Description 

Debug option for additional print 
default 0. Off, or 1. 

Number of blades 

Propeller diameter, feet 

Inner-most blade station at 50% chord 

Defined in input as SWPOPT 

Defined in input as THTDES 

Defined in input as THTCUT 

Defined in input as ZKCUT 

Station radius/blade radius 

Maximum blade thickness/blade chord 

Blade chord/propeller diameter 

Design lift coefficient 

Blade twist, from plane of rotation, 
degrees 

Defined in input 
Defined in input 
Defined in input 
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AIRFOIL TYPE 


Defined in input as BAFL 
NO. OF % CHORD Defined in input as ZNPCOV 

LIST OF % CHORD Defined in input as PCTCHD 


NO. OUT STATIONS Defined in input as ZNIS 

LIST - STATIONS Defined in input as ZBLDST 

AERODYNAMIC OUTPUT 

AEROEXEC : Figure 22A summarizes the flight condition, blade 

geometry and other parameters that are used in the panel code 
performance calculation (AEROEXEC) . Definition of the parameters 
in Figure 22A follows: 

Parameter Description 


TEMP, DEGF 
RHO/RHO STD 
SPEED OF SOUND 
ADVANCE RATIO 
FLIGHT MACH NO. 
FLIGHT SPD FPS 
RPM 

TIP HEL. MACH 

START BLENDING 
DIAMETER 
NO. BLADES 
NO. INPT STA. 
FREQ. OF UNST. 
NO. NODAL DIA. 

K - DOWN 


Ambient temperature, °F 
Ambient density/standard density 
Ambient speed of sound, fps 
Flight velocity/ ( rps*diameter ) 
Flight Mach number 
Flight velocity, fps 
Propeller, rpm 

Mach number based on resultant 
flight and tip speed velocities 

Defined in input as QMBLEN 

Propeller diameter, feet 

Blade number 

Number of input stations 
Defined in input as QQ 
Defined in input as QK 
Defined in input as QKDOWN 
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K - START 
INPUT STATIONS 
B/D 


Defined in input as QKSTART 
Blade radial stations/blade radius 
Blade chord/propeller diameter 


TOTAL TWIST 
MCA/D 


FA/D 


T/B 

CLD 

SWEEP 


Operating twist from plane of rotation, 
degrees 

Distances from pitch change axis 
to blade mid-chord point along the 
helix/propeller diameter (ref. 

Fig 3.1 of Vol.I) 

Perpendicular distance from 
helix to the intersection of 
the blade mid-chord and 
mid-camber point/diameter (ref. 

Fig 3.1 of Vol.I) 

Max. blade thicknesses/blade chord 

Blade design lift coefficients 

Sweep angles between 50% mid chord line 
and resultant inflow velocity 


Another output table is given in Figures 22B and C. This is a 
table of the camber height/blade chord as a function of local 
radius/blade radius and local chord/blade chord which has been 
calculated from the input in BLADEGEO. The propeller driveshaft 
power which has been converted to acoustic energy has been 
calculated and is listed on figure 23 in coefficient form with 
the heading "SOUND POWER". 

LSTBUG = 1. also prints intermediate non-linear iteration steps, 
as shown in Figures 23. The iteration continues (and has been 
omitted herein) until CPFRONT, the power coefficient ceases to 
change. Figure 24. The parameters shown are: 


SOUND POWER = acoustic power loss/ (density* (rev/ unit time) 3 * 

(diam) 5 ) 

CP = power / (density* (rev/unit time) 3 * (diameter) 5 ) 

CT = thrust/ (density* (rev/unit time) 2 * (diameter) 4 ) 

EFFICIENCY = apparent efficiency = CT*advance ratio/CP 
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VORTEX OUTPUT 


VORTEXEC: The output from the vortex performance calculation is 
given in Figure 25 and includes four pages of output described 
below. 


The first page, Figure 25A defines the selected operating 
condition: 


Parameter Description 

FLIGHT The flight Mach number 

MACH # 


TIP ROT. The tip rotational Mach number 

MACH # 


ADVANCE The advance ratio 

RATIO 


TEMP. Ambient static temperature, °F 

DEG F 


DENSITY 

RATIO 


Density ratio, ambient density/standard 
level density 


# RADIAL The number of output radial stations 

STATIONS 


# CHORD- The number of chordwise stations, QNCP 

WISE 
STATIONS 


# SPAN- The number of spanwise modes or control 

WISE points, QSMN 

MODES 


MODOPT The spanwise mode type option, defined in 

LSTPARMS input as QMODOP 

NBOPT The lead-edge panel load shape option, 

defined in LSTPARMS input as QNBOPT 

The second output page Figure 25B describes blade geometry 
parameters and some force coefficients from the potential flow 
lifting surface calculation. A brief description of the 
tabulated parameters follows. 


Parameter 


Description 


THETA 3/4 


The blade angle at the .75 radius 
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# Blades 
DIAMETER 
SPINNER CUTOFF 

RADIAL STATION 

CHORD/ DIAMETER 

THICK/ CHORD 

MID-CHORD COORDINATES 

XMC 

YMC 

ZMC 

FACE ALINE 

MID-CHORD ALINE 

MID-CHORD SWEEP 
L.E. SWEEP 
ALPHA 2-D 
ALPHA 3-D 

L.E. CAMBER 


Number of propeller blades 

Propeller diameter, feet 

Innermost airfoil shaped blade 
radius/propeller radius. 

The selected output spanwise radial 
station/propeller radius. 

Local chord/propeller diameter for each 
radial station. 

Local section maximum blade thickness 
local blade chord. 

The coordinates, at 50% chord are defined 
in figures 5, 6 & 7. 

Local radial coordinate/propeller radius. 

Local inplane coordinate/propeller radius. 

Local axial coordinate/propeller radius, 
(positive downstream) 

Perpendicular distance from the local 
mid-chord location to the local helical 
plane/propeller diameter. 

Distance along the local helical plane to 
the local mid-chord location/propeller 
diameter. 

Local sweep angle between the mid chord 
line and the resultant velocity 

Local sweep angle between the lead edge 
and the resultant velocity, degrees. 

Two-dimensional angle of attack used to 
determine airfoil tabulated C 0 , degrees. 

Angle difference between the twist angle 
and the advance angle, both defined below, 
degrees . 

Lead edge camber angle (angle between 
camber line at the lead edge and the chord 
line) , degrees 
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L.E. A.O.A Lead edge angle of attack 

(ALPHA 3-D minus camber angle), degrees. 

ADVANCE ANGLE Angle between the resultant velocity 

neglecting induced flow and the local 
rotational velocity, degrees. 

TWIST ANGLE Local blade angle measured from plane of 

rotation, degrees. 

Note : all lift and drag coefficients are 
non-dimensionalized by qr*b*dr where : 


CDLST 


CDMIN 


DO 

CDCLVD 


qr - dynamic pressure of the resultant 
velocity 

b - blade chord defined above 

dr - unit spanwise length 

The drag coefficient obtained from the 
airfoil tables at ALPHA 2-D 
Where drag is a downstream force along the 
advance angle 

Minimum drag coefficient. 

Drag coefficient at lift coefficient = 0 

Drag to lift ratio 


CDUAP2 


Total drag coefficient equal to : 
CDLST + CDCLVD * CLLST 


CLLST 


Lift coefficient calculated from lifting 
surface 


L.E.K., MAG Flat plate load coefficient where, 

L.E.K., MAG = del Cp * (fc/(l-fc) )**.5 
del Cp = the differential pressure 
coefficient between blade upper and lower 
surfaces and 

fc = fractional chord length 

L.E. THRUST Component of leading edge suction parallel 

to leading edge camber line equal to 
Thrust/qr*b*dr 

qr - dynamic pressure of the resultant 
velocity 

b - blade chord defined above 

dr - unit spanwise length 
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RELATIVE MACH # 


Local advance Mach number = 
(M 0 **2+M T **2) **.5 
Mo - freestream Mach number 
Mt - tip speed/speed of sound 

DESIGN CL Local blade design lift coefficient 

(camber) 


The next page, Figure 25C, of output present the lift and drag 
coefficient components associated with the vortex generated lift, 
as a function of local radial station. The lift force is 
perpendicular to the local chord; the drag force is parallel to 
the local chord, (as described in Volume III appendix) . The 
components are : 


Parameter 
POT LIFT 

L.E. VORTEX 
TIP VORTEX 

AUGUMENTED VORTEX 
TOTAL LIFT 


Description 

Calculated lift coefficient from lifting 
surface theory, (same as CLLST) 

Normal force coefficient from the lead 
edge vortex which occurs if L.E. sweep is 
greater than 12 Deg. 

Normal force coefficient which occurs 
because of flow separation in the blade 
tip region causing the tip vortex to 
reform on the blade tip upper surface. 

Normal force coefficient from vortices 
shed from the lead edge vortex passing 
over the blade tip. 

Total lift coefficient 


POT + LE VORT 


POT + SE VORT 


POT + AUGVORT 


POT DRAG 


Sum of the calculated lifting surface lift 
coefficient and the lift coefficient from 
the lead edge vortex. 

Sum of the calculated lifting surface lift 
coefficient and the lift coefficient from 
the tip vortex. 

Sum of the calculated lifting surface lift 
coefficient and the lift coefficient from 
the augumented lift. 

Calculated drag coefficient from lifting 
surface theory ("induced" or "vortex" 
drag) 
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POT + LE VORT 


Sum of POT DRAG and drag coefficient from 
the lead edge vortex 

POT + SE VORT Sum of POT DRAG and drag coefficient from 

the tip vortex 

POT + AUG VORT Sum of POT DRAG and drag coefficient from 

augmented normal coefficient 

TOTAL DRAG Total drag coefficient 

The final page. Figure 25D, of AEROEXEC output shows the spanwise 
change in power coefficient (DCP/DX) , thrust coefficient (DCT/DX) 
and efficiency (DETA/DX) for each of the components in lift and 
drag discussed above, and then finally the integrated, power 
coefficients, thrust coefficients and efficiency. 

WAKE OUTPUT 


WAKEEXEC ; The output from the wake velocity calculation Figure 26A 
presents a table of the following velocity parameters, defined in 
the velocity diagram of Figure 26B. 


Parameter 

PHI 

Vo 

Uo 

U 

W 

Vx 

VT 

RADIUS 

XMBAR 

POTENTIAL 


Description 

Cylindrical coordinate angle, degrees. 

Zero degrees at the sheet. 

Freestream velocity 

Local resultant velocity. 

Induced velocity in resultant velocity 
direction (+ direction shown in diagram) . 

Induced velocity perpendicular to 
resultant velocity direction 
(+ direction shown in diagram) 

Induced velocity in axial direction 
(+ direction shown in diagram) . 

Induced velocity in tangential direction 
(+ direction shown in diagram) . 

Local radius/propeller radius at which the 
wake velocity is calculated. 

Axial distance/propeller radius from the 
propeller pitch change axis (+ downstream) 

Potential flow velocity component. 
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VISCOUS 


Viscous wake velocity component, for this 
case the viscous wake is very thin, and 
appears only at PHI =0. The viscous wake 
component is centered at PHI=0, is thin at 
the airfoil trailing edge, and the width 
increases with downstream distance. 

TOTAL Sum of the potential and viscous velocity 

components including freestream velocity. 

NOISE OUTPUT 


N0I2EXEC : The first page of the noise program output (see Figure 
27A) summarizes the operating conditions and observer positions 
selected by the user. The first few lines repeat the header 
information. The operating conditions are then listed : 


Parameter 

FLIGHT MACH NUMBER 
TIP ROTATIONAL MACH 
TIP MACH NUMBER 

TEMPERATURE, DEG F 
SPEED OF SOUND 
HZ BPF 

PROPELLER DIAMETER 
NUMBER OF BLADES 
PAMB/2116 

RPM 


Description 

Flight mach number 

Tip speed/speed of sound 

Resultant of flight and tip rotational 
Mach numbers 

Ambient temperature, deg.F. 

Ambient speed of sound, fps 

Blade passing frequency. Hertz (cycles 
per second) 

Propeller diameter, ft. 

Number of blades on propeller 

Local atmospheric pressure/sea level 
standard pressure 

Propeller speed, revolutions per 
minute 


The calculation instructions list some of the parameters 
important to the noise calculation. They are: 


Parameter 

ROOT RADIUS RATIO 

TIP RADIUS RATIO 


Description 

Most inboard location on blade 

Furthest outboard location on blade 
(normally 1.0) 
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NUMBER OF HARMONICS - Maximum noise harmonic order to calculate 

(Expression) PTS CHORDWISE - Number of chordwise integration 

mesh points 

DIRECTIVITY POINTS ANALYZED - Number of axial directivities 
NUMBER OF RADIAL STATIONS ON BLADE 

(Expression) PTS/CYCLE RADIAL INT - Number of radial integration 

mesh points 

(Expression) PTS/CYCLE W INT - Number of frequency integration 

mesh points in the range 
1/1+MX to 1/1-MX (used in the 
near field calc, only) . MX is 
the flight Mach number. 

W EXTEND POINTS - Number of points in the frequency 

integration outside 

the range 1/1+MX to 1/1-MX 

NEAR FIELD ANALYSIS INDICATOR - 0 for far-field, 

1 for near-field calc. 

DETAIL PRINT INDICATOR - Number of harmonics for which a detailed 

print will be obtained (set by ZIPDET) 

QUADRUPOLE TERM CALC. CODE - Not 0 to calculate "instant" 

quadrupole noise. 

The observer directivity points lie on a line parallel to the 
propeller rotation axis. The distance between the axis of 
rotation and the axial observer line is the altitude or sideline 
distance. The relationships between the visual and retarded 
distances and angles are shown in Figure 27B. 

The addition of the boundary layer and wake displacement 
thickness to the airfoil thickness is indicated by the next line. 


AZIMUTHAL DIRECTIVITY ANGLE - see Figure 27C for a description of 

this variable. 

On the next page (Figure 28) is a header indicating the beginning 
of the noise calculation results. 
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At this point a description of the detailed printout (obtained 
when ZIPDET is not 0) is appropriate (see figure 29A) . For each 
harmonic up to ZIPDET, a table of information is presented as a 
function of radial position. This includes the following 
information: 


Parameter 

RADIAL STATION INDEX 

RANGE 

KX 

KY 

PHASE LAG DUE TO 


Description 

Counter for radial elements 

Inner and outer limits of radial element 

Chordwise wave number 

Normal -to-chord wave number 

Blade sweep 


PHASE LAG DUE TO OFFSET 


SUM OF SOURCES 

MONOPOLE 

DIPOLE 

QUADRUPOLE 

AMP 

PHASE 

RADIAL LOAD 


From the helical surface 

Total noise including radial load and 
"instant" quadrupole 

Blade thickness (plus boundary layer and 
wake displacement) sources. 

Lift, drag, and radial loading sources. 

Nonlinear sources ("instant" quadrupole) 

Linear noise amplitude scaled to an 
arbitrary value 

Relative phase of noise components 
Noise due to radial forces at blade tip 


BLADE TOTALS PASCALS (RMS) Sound pressure levels 

DB - NO INTERFERENCE Noise level neglecting the effects of 

sweep and offset. This is computed to show 
how much noise reduction is due to blade 
sweep and offset. The noise reduction 
depends on the operating condition as well 
as the propeller configuration. The noise 
level neglecting blade sweep and offset is 
computed by integrating the noise 
contribution from the spanwise strips 
using the phase angles PHIO and PHIS 
(offset and sweep) . 
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DECIBELS Resultant noise level in dB, including 

sweep and offset. 

The following page (Figure 29B) is a summary of the noise 
calculation results for the first directivity. This type of page 
is repeated for as many directivities as were requested by the 
user. The information presented is: 


X 


Retarded axial position 


XI 


Visual axial position 


THETA 

HARMONIC ORDER 
FREQ HZ 

SUM OF 3 SOURCES 


DB 

PASCALS 

PHASE 

MONOPOLE 

DIPOLE 

QUADRUPOLE 

RADIAL LOAD 
INTEGRATION MESHES 
JAMX 
KK 

KEXTND 


Retarded observer angle 

Multiple of blade passage frequency 

Harmonic frequency, Hertz 

Indicates the total noise for all sources 
(including blade thickness, boundary layer 
and wake, lift and drag, quadrupole 
and radial blade loads) . 

Noise level. Decibels 

Sound pressure level in Pascals 

Relative phase of the noise component 

Noise due to blade thickness plus 
displacement thickness 

Noise due to lift and drag 

Noise due to nonlinear sources (normally 
not used herein) 

Noise due to tip edge radial forces 

Number of points used in integrations 

Number of radial mesh points 

Number of frequency integration points in 
range from 1/1+MX to 1/1-MX, where MX is 
the flight Mach number. 

Number of points outside the KK range 
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The FINAL NOISE SUMMARY TABLE (Figure 30) is presented on the 
next page. The header is repeated at the top of this page, 
followed by the sideline distance and a table of the axial 
observer locations. Then the total noise and its components are 
listed for each directivity. 
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SECTION 5 

SYSTEM DESCRIPTION 
Language Req uirements 

The source programs included in the delivered code are written in 
the FORTRAN 77 language. This code was originally written on the 
IBM 3090, and compiled using the IBM LEVEL 1.4.1 (May, 1985) VS 
FORTRAN compiler. Extensive use of double precision variable 
typing in the form of DOUBLE PRECISION, REAL * 16, IMPLICIT 
DOUBLE PRECISION and COMPLEX * 16, along with the implied double 
precision arithmetic was reguired on the IBM 3090 to achieve the 
numerical precision demanded by the analysis methods. The 
delivered code contains these type statements. However, on the 
CRAY these double precision statements are "turned off" by use of 
a compiler option. 

Note: REAL *16 and COMPLEX * 16 are not Fortran 77 

Statements . 

Library Requirements 

The IMSL (International Mathematical and Statistical Library 
Version 10) , libraries MATH and SFUN are reguired by the code. 

The IMSL single precision library is used by the CRAY version, 
while, the IBM 3090 requires the IMSL double precision libraries. 

The standard FORTRAN library functions are required by this 
code, and have not been delivered with the code, e.g. COS, SIN, 
etc. 

Subroutine CLOCK, and DATE are CRAY calls to obtain system time 
and date respectively. 


File Allocations 

The file requirements are shown in figure 31, "UAAP System File 
Requirements". Conceptually these files are used in the 
following manner: 

5 Input options, specifications, and commands are 
entered on file 5 by the user. 

6 Program output file for hardcopy printing. Contents 
of this file are detailed in the "OUTPUT 
Description" section of this manual. 

8 The major amount of CPU time spent by the program 

involves generating and inverting a matrix, K* 1 , in 
the lifting surface aerodynamic section of the code, 
LSTALS. The K' 1 matrix is written to file 8, and 
later re-read by the program. If the user wishes to 
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save this file after program termination, the user 
will have to create appropriate machine dependent 
Job Control Language code to do so. Saving of the K' 1 
matrix will allow other cases to be run with this K' 1 
thus saving CPU time. 

9 This file contains the "THICKNESS VECTOR" generated 

by the lifting surface aerodynamic section of the 
code, LSTALS, it will be read into the program 
later. As with the K’ 1 matrix on file 8 the user may 
save the "THICKNESS VECTOR" for future use. 

11 Subroutine INPTRW reads the entire file 5 and echo's 

it to the output file 6, and to a temporary storage 
file 11. File 11 will be rewound and re-read by 
individual subcommands. 

44 This file is used to transfer data between the Panel 

portion of the aerodynamic code and the Vortex 
calculation portion of the aerodynamic code. 

50 This file is used to transfer data between the 
aerodynamic and the noise sections of the code. 

51 Scratch file, for use by the acoustic portion of the 
code. 

52 Scratch file, for use by the acoustic portion of 
the code. 

Input Data Storage Flow Chart 

Input data is arranged according to the functional area in which 
it will be used. Each group of input is read into the functional 
area in which it is to be used, and is stored in a Fortran 
Labeled Common. Figure 32 illustrates the functional areas into 
which input is read, and the name of the Labeled Common in which 
the data is stored. With the exception of BLDGEM, and VELGRD, no 
substantial amount of calculations are performed on the input 
data before storage in the Labeled Common areas. In BLDGEM the 
input blade geometric description is transformed into a form 
useable by the other sections of the code, before storage into 
BLDGEO Labeled Common area. In VELGRD, the input axissymmetric 
velocity gradients are normalized by the freestream velocity 
before storage in the VII 1C1 Labeled Common area. 
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TREE Structure 


A tree structure which identifies the calling sequence in the 
program has been created. The tree structure provides the 
calling sequences possible during execution, but program options 
which decide the actual calling sequence during execution are not 
accounted for. 

The number on the left of the tree. Figure 33, indicates the 
"level" at which the subroutine (whose name appears next to the 
level nimber) is called. A brief description of the major 
subroutines is provided adjacent to the subroutine name. 

For example, referring to FIG 33, the 

MAIN program calls DATE to get the current date, this is a "level 
1" call : 

MAIN calls CI/DCK to get the current time of day, 

MAIN calls INPTRW to read file 5 and echo the input to 
file li. 

INPTRW calls NEWPGI to print a page heading, this is a 
"level 2" call, 

MAIN calls BLDGEM which is the top level routine for 
calculating the blade geometry. 

BLDGEM calls EXIT, a System Routine. This call occurs only 
when a particular error occurs. 


Subroutine References 

Each subroutine contained in the delivered code has been 
alphabetically listed in figure 34. Adjacent to the subroutine 
name is brief description of that subroutine if one exists in the 
source code. All of the subroutines referenced by that 
subroutine are listed on the next line adjacent to the word 
"CALL(S) :". 

Figure 35 provides an alphabetized list of subroutines, followed 
by a list of subroutines which call them. 

List of Su broutines 

Figure 36 provides an alphabetized list of subroutines which 
include labeled common areas, and the names of the labeled common 
areas. Subroutines names appear on the left while labeled common 
areas appear to the right of the colon. Figure 37 is a cross 
reference of labeled common blocks and subroutines using them. 
Common block names appear on the left while subroutines in which 
the common block is used appear to the right of the colon. 
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SECTION 6 
ERROR CODES 

This section contains a description of the printed error 
messages, along with possible causes and remedies in alphabetized 
order . 

"AIRFOIL NO. = ... OFF AIRFOIL DATA AND IOFF ..." 

The program obtains profile drag from airfoil data 
tables as a function of lift, mach number, airfoil 
thickness ratio, and airfoil chord. This message 
indicates one or more of the above parameters 
exceeded the limits of the stored data. 

Recommendation : check the output for reasonable 
values of drag coefficients. Use AIRPARMS inputs 
CDMULT and DCD to correct unreasonable values. 

"AIRFOIL . . . OFF AIRFOIL DATA " 

Same as above. 

"BDSO ... " 

There are a number of errors, generally resulting in 
program termination, which are associated with the 
blade geometry generator; these messages are 
prefaced with BDSO. They arise because the blade 
surface generated from user input in BLADEGEO was 
not sufficient to allow intersection with the output 
stations (conical surfaces) specified in BLADEGEO 
variable ZBLDST. Solution; The problem can usually 
be avoided (with only a slight loss in program 
accuracy) by setting THTDES = THTCUT in BLADEGEO, 
and assuring that X(l) < ZBLDST (1) and X(10) > 

ZBLDST (ZNIS) in BLADEGEO input. 

••*** CONTROL POINT STATIONS NOT WITHIN ..." 

The radial difference between control points, QZAR 
in LSTPARMS , and ZBLDST in BLADEGEO > . 02 . 

Solution : Change the either QZAR or ZBLDST input 
values. 

"FAILED TO CONVERGE ON CL ..." 

An iteration failed when trying to obtain the 
profile drag data from the airfoil tables. 
Recommendation : check the output for reasonable 
values of drag coefficients. Use AIRPARMS inputs 
CDMULT and DCD to correct unreasonable values. 

"FOR OUTPUT POINT NUMBER ..." 

Interpolation of the flowfield from VELGRADS is 
performed in two ways ; radial first, axial second, 
and then axial first, radial second. Usually these 
interpolations yield the same value of velocity at a 
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fixed radial and axial location. However, for this 
point the difference in interpolations exceeded 1 %. 
Probable cause : the flowfield is not smooth in 
either or both the radial and axial directions. 

"FOR OUTPUT POINT NUMBER ... THE AXIAL ..." 

The axial extent of coordinates input into VELGRADS 
was not sufficient to allow interpolation of the 
flowfield at either the blade leading or trailing 
edge. 

Solution : add more points to the flowfield input in 
VELGRADS in the axial direction. 

"FOR OUTPUT POINT NUMBER . . . THE RADIAL ..." 

The radial extent of coordinates input into VELGRADS 
was not sufficient to allow interpolation of the 
flowfield at either the blade root or tip. 

Solution : add more points to the flowfield input in 
VELGRADS in the radial direction. 

"INVALID INPUT CHARACTER ..." 

The load routine has found a character in column 1 
which is not an "L", "C", "E", or blank. 

Possible causes : "END" record omitted; numeric data 
in column 1. 

"THE INPUT CONTAINS AN UNKNOWN COMMAND ... " 

A COMMAND was expected at this point in the input, 
however, the 8 printed characters do not represent a 
recognizable commnand. 

Possible causes ; mis-spelled command; no "END" 
record in load format. 

COMMAND didn't start in column 1. 

"THE INPUT CONTAINS AN UNKNOWN SUBCOMMAND.." 

A subcommand was found which didn't match the 
acceptable subcommands. 

Possible causes ; incorrect spelling; incorrect 
placement of parenthesis. 

"THE NUMBER OF HEADER CARDS " 

More than 10 header records were encountered during 
processing of the HEADER command. 

Possible causes ; more the 10 records following the 
HEADER command; no "END" record after the last 
header record. 

"V1S1X1 DIMENSIONED 20X20 ..." 

The flowfield input in VELGRADS was too large; the 
maximum number of coordinates in the axial or radial 
direction is 20. 
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Solution : reduce the number VINRD and/or VINAS to 

20 . 

****** MM2 AND MM3 MUST BE LESS THAN 100 ..." 

See QMM2 and QMM3 limits in LSTPARMS. 

"***** NBOPT MUST BE 1 OR 2 ..." 

See QNBOPT limits in LSTPARMS. 

"***** NCP GREATER THAN MAX NCP ..." 

See QNCN limits in LSTPARMS. 

"***** NIS GREATER THAN MAX NIS ... " 

See ZNIS limits in BLADEGEO. 

"** NOIZEXEC CALLED WITH KEYWORD ..." 

An unrecognizable subcommand was found. The 
subcommand must be "EXECCASE" . 

****** NSM GREATER THAN MAX NSM ... " 

See QNSM limits in LSTPARMS. 

"***** NSM * NCP GREATER THAN ..." 

Ensure that QNSM * QNCP < 1000 in LSTPARMS. 

"*** NUMBER OF NODAL DIAMETERS (K) MUST ..." 

If QQ = 0 then QK must equal zero in the LSTPARMS 
input. 

"***** NX GREATER THAN ..." 

See ZNPCOV limits in BLADEGEO. 

••**** THE DATA FOR THIS RUN & THE DATA ..." 

A label is attached to the "K-INVERSE" matrix 
describing the parameters used to create it. The 
current value of the input variable listed below the 
error message does not agree with that in the 
attached label. 

Cause : This generally arises because a "K-INVERSE" 
matrix from a previous run was used as input to this 
run (see QPART1 in LSTPARMS) and the data used to 
generate that matrix is different than that being 
used to run the current case. 

"*** THE NO. OF DIRECTIVITY POINTS ..." 

See ZNX limits in NOIZPARM. 

"*** THE NO. OF HARMONICS IS GREATER ..." 

See ZMMAX limits in NOIZPARM. 

"*** THE NO. OF HARMONICS IS GREATER THAN ALLOWED 150' 

••*** THE NO. OF HARMONICS X THE NO. OF BLADES IS LIMITED TO 
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1000 ' 

"*** THE NO. OF RADIAL INTEGRATION POINTS IS LIMITED TO 2000 
"*** THE NO. OF DIRECTIVITY POINTS IS LIMITED TO 20' 

Cause - input values larger than permissible 
Solution - correct input 

" KKR2 VALUE = computed LARGER THAN HANKS DIMENSION FOR HARMONIC 
number, MMAX IS BEING RESET TO PREVENT ERROR ' 

Cause : limited storage in Hankel function array and 

conservative selection of integration step size. 

This error is most likely to occur at high flight 
Mach number , forward directivities, and/or high 
harmonic order. 

Solution : 

1) Reduce distance forward of plane of rotation. 

2) See note on selection of ZKXTND in the input 
section. 

3) Reduce ZKMIN and ZOMEGA from the default values 
of 10. Values of 6 or 7 can give acceptable 
results. 
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Input data 
(steady) 
(calculation) 


Input data 
(unsteady) 
(calculation) 



UAAP system overview 
Figure 1 
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Figure 4 Blade geometry nomenclature (2-D coords) 
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Figure 5 Blade geometry nomenclature (RXY coord) 


MFA?1?4 DO 4 


Max. thickness (THKOB 0 CHORD) 



Figure 6 Blade geometry nomenclature (RXY coord) 



Figure 7 Blade geometry nomenclature (RXY coord) 


A Z (Pitch change axis) 



Figure 8 Blade geometry nomenclature (XYZ coords) 


READ / PRINT ALL INPUT DATA i INPTHR INPlfTRH TIME :HH:MM:S$ DATE : MM/DO/YY UAAP — NASA 






O O O' O CO pH 

XI O CM >6 pH HI 

XI O U| rH »H sO 


rH XI O 

CM O' X* 

^ n© H 


CM >A CO H sD 

H H 9 H H 


cm 9 - cp h 9 1 

H H H H 


O nP O' 

rH pH 9 1 

I CM 


O O CO 

2 S S 


O O O' 

N U) M 

9* CM XI 


9" U| U| 

fh in i 

M 


CM O' © CM 

fH 9’ rH fH 


* 

t 

o 

4 

CM 

o 

O' 

i H 

O C 

i © 

d 

O 

d 

o 

d 

d 

rH 

d 

K 

pH 

d 

4 

pH 


d 

pH 

XI 

in 

rH 

-1 

1 

rH 

9 

CM 

rH 

1 

rH 

9 

CM 

H 

1 

pH b 

C 

3 7 

pH 

3 

t 

pH 

9* 

CM 

i 

pH 

in 

CM 

1 

pH 

3 

1 

rH 

in 

CM 

1 

pH 

in 

CM 

1 

pH 


O' O O O 
■9 co o o 
M N N O 

rH 9* d 9 
M Kl vt H 
CM CM I 

IflOOO 
>0040 
H U) 9" O 


•H O O O 
NUiroo 
O N CM O 

9 CM CO 9 
CM CM 9 H 
CM CM I 

XI O © O 
9 N S O 
H Ifl (O O 

H M N 9 
N H 9 H 
CM CM I 

CM O O O 
N M 9 O 
in CM CM o 


O (O 4 O I 

O N 4 4 I 

ONMWl 

9 «0 H 4 I 

I H IO 9 « 

I CM 1 

0(0901 

O 9 in 4 I 

O rt H N I 

9 N N 4 ' 

I rH CM >9 

I CM < 

>091001 

>otomo 

i O CO O CSJ I 


© CM CM O 
0900 
Ifl O 9 9 

l CM 9 rH O 

i h xi in 


O O O <9 I 
00004' 
CM © O pH 

O H H 9 i 
Ifl H H I 
CM I I 


I O O © *H 
i o m o o' 
■ sO CM in CM 

10099 

\ m h i 

I CM I 


rH © O CO 
CO H 4 O 
9 9 fh m 

H H H 9 

xi in in > 

CM CM 


OC0SOOO99inOOp4f0OOOin9N4OOM9C0OOOON4 

9e0OC0O99KINHin«9NKI(0OH9HinH4999NONS9N 

LhKunH94NinincouiN44in9N99Mn9a)OKiNinin9H4«o 

N9HOH(049H9HMn9HOH499HOH499HKIHinKI9 

i m m in i i xi m m i i xi m> m i i xiuSiniT x> in in i i 

CM CM CM CM CM CM CM CM CM CM 


omoooooHioooooHinooioo>o9ooinxiH>o 

ONOMMOO(09H9inOMHK99999©MOO9N 

:in4CJN(OOON(M9^Ntn9MKI99mOrONKUnKIHM 


CM CM CM CM CM CM CM CM 

oom90HooinmM9oo9«HNooin9o 

KH4ine09H»nNOC0HinN©9«0MN499(n 

MN44MMWffl(0OW9NO9M99inOH4in 


N 9 9 H 

CM in 4* fH 


9 m 9 »^ > 
rH CM 9 


O O O' o < 
O O XI CM < 
O © CM CM 4 

99HM 
fh cm 9 ' 

I CM 1 

o o m o i 
a xi xi 9 < 

O H 4 N ' 

9 9 CO (0 I 
l rt MO’ 

I CM l 


m o o o < 

H 9 XI O < 

H 4 9 O < 

9 O 9 H t 

N 4 9 H i 

CM CM I 

r*- o © o < 

HMOOI 

i xi in o O 4 

I H (0 9 H i 

CM Ml 9* M * 

CM CM I 

>90004 
IM990I 
i N 9 4 O l 

i co in «o m < 

H N 9 H i 

CM CM I 


HNN4009NM009MONOHa4(ON 

h cm in m m i cm in m i i cm in m i » cm 

I CM CM I CM CM CM CM 

OH9OOO>0MinOOM9inOOQN9MH 

O4in49in9^9mHNHinCDinH9MC0N 

O4HKIONinC0N4NM49W9HWMOXI 


CM © © O' 

i so in m i 

CM CM I 

• 0009- 

i h 4 in SI 

I 00 N CM 4 

19 9 0 9 
I49H I 
CM CM I 

' © © O XI 
|N CM m H 
. 4 9 N 9 

i in O' © o 
I N 9 f4 • 
CM CM I 


49WO9N499HC04N9 

cm m m i i cm sf in i i <m 

CM CM CM CM 

S90040NUn00NHH(D 

O4O9OHNH9HNC0Hin 

O9HMNOHint0«0HNNin 

4HCpO9(04H00O(D44H 

cm in m i i n 9 in i i cm 

CM CM CM CM 

ooooocookooopmo9 

in9fflXI9lfl494MinNOQ 

m«99 H 9 49NWOKINO 


HHMnC0KNHsDlfl(0N9CMlfl9«NMMinWC0 

4 in i i cm vo m i i cm in m i i cm in in i i 

CM CM CM CM CM CM CM CM 

OOW99KIOOC0N9NOOOf-(49OOOlflO 

O9NS49NM9H(0NNO94H9tMNNOl0 

90in4HMHN40NKlK9®XI9M9N0Mn 

FMCMhI|AK<^CMCMsOinr^>^C0CMin9 I K<^KICMinNlK 

4 io i I cm 4 in i i cm in m i > cMinmii 

CM CM CM CM CM CM CM CM 

OO99N«OOO9N9OOrtC0NXIOO«99 

OH4940<p9(04HNMOU|4NOIA4NN9<0 

MU1MONUiinMAM94O9NlA44999(0N 

9NN44H9N4ui4HKNIA94H9N9KI4 

m m i t csimmii cm w m i i cm in m • i 


i m m i t cm in m i 

CM CM CM CM 


cm m m i 

CM CM 


m«999(0lfl(0NO «Nlfl9 
H II) 9 I I PM 9* m I I rH 
CM CM CM CM 


ooaO'ocM©oof-400'ooom©soooo©o 

HlANC099NM9H>0(0NOHNirMM9l0Kl9«0 

U)HNU|NON990909U1N0HH9lfl9Nin 

^^N4lfl94^4uiwi4^9fOinin49lflN994 

in in t i h in in i i nmin i i ph m m » i 

CM CM CM CM CM CM CM CM 


Figure 10 Program Output Listing : Input Echo 
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Figure 11 Program Output Listing : Input Echo 
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Figure 12 Program Output Listing : Input Echo 
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Figure 13 Program Output Listing : Input Echo 
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Figure 14 Program Output Listing : Runparms Output 
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Figure 15 Program Output Listing : Noizparm Output 
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Figure 16 Program Output Listing : Vortparm Output 
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Figure 17 Program Output Listing : Lstparms Output 
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Figure 18 Program Output Listing : Airparms Output 


PROPELLER HAKE ANALYSIS INPUT :HAKPRM HAKEPARM TIME :HH:MM:SS DATE ; KM/DD/YY UAAP — NASA 
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Figure 19 Program Output Listing : Wakeparm Output 
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Figure 20B Program Output Listing : Velgrads Output (continued) 
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Figure 20C Program Output Listing : Velgrads Output (continued) 
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Figure 21 Program Output Listing : Bladegeo Output 
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Figure 22A Program Output Listing : Aeroexec Output 
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Figure 22B Program Output Listing : Aeroexec Output (continued) 
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Program Output Listing : Aeroexec Output (continued) 
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Figure 23 Program Output Listing : Aeroexec Output-Performance 




Figure 24 Program Output Listing : Aeroexec Output-Performance 
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Figure 25A Program Output Listing : Vortexec Output 
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Figure 25B Program Output Listing : Vortexec Output (continued) 
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Figure 25C Program Output Listing : Vortexec Output (continued) 
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Figure 25D Program Output Listing : Vortexec Output (continued) 


WITH CORRECTION FOR NON-LINEARITY 
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Figure 26A Program Output Listing : Wakeexec Output 


Figure 26B Wake exec output description 
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EXECUTE NOISE CALCULATION : NOZCLC NOIZEXEC EXECCASE TIME :HH:M 1 :SS DATE : MM/DD/YY UAAP — NASA 
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Figure 27B Visual and Retarded Location Relations 
(x and are positive in this example) 




Unsteady loads are defined with respect to the reference blade 
position at time T = 0. The observer azimuthal position is defined by 
the angle ^measured in the direction of rotation from the reference 
blade position. 


Figure 27C Azimuthal observer angle. 
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Program Output Listing : Noizexec Output 


BLADE ELEMENT RESULTS FOR HARMONIC NUMBER 1 
FREQUENCY IN HERTZ 2 617.4, DIRECTIVITY POINT 
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Figure 29A Program Output Listing : Noizexec Output 


OUTPUT DIRECTIVITY RESULTS : NOZSMS NOIZEXEC EXECCASE TIME :HH:MM:SS DATE : MH/OD/YY UAAP NASA 




Figure 29B Program Output Listing : Noizexec Output (continued) 


BLADE ELEMENT RESULTS FOR HARMONIC NUMBER I 
FREQUENCY IN HERTZ* 617.4, DIRECTIVITY POINT 
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Figure 29C Program Output Listing : Noizexec Output (continued) 


OUTPUT DIRECTIVITY RESULTS : NOZSMS NOIZEXEC EXECCASE TIME :HH:MH:SS DATE i MM/OO/YY UAAP NASA 
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Figure 29D Program Output Listing : Noizexec Output (continued) 


BLADE ELEMENT RESULTS FOR HARMONIC NUMBER 1 
FREQUENCY IN HERTZ® 617.4, DIRECTIVITY POINT 
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Figure 29E Program Output Listing : Noizexec Output (continued) 


OUTPUT DIRECTIVITY RESULTS i NOZSMS NOIZEXEC EXECCASE TIHE :HH:hW:SS DATE : MM/OO/YY UAAP — NASA 
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Figure 29F Program Output Listing : Noizexec Output (continued) 


BLADE ELEMENT RESULTS FOR HARMONIC NUMBER 1 
FREQUENCY IN HERTZ* 617.4, DIRECTIVITY POINT 
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Figure 29G Program Output Listing : Noizexec Output (continued) 


OUTPUT DIRECTIVITY RESULTS : NOZSMS NOIZEXEC EXECCASE TIME :HH:MM;SS DATE : TW/OD/YY UAAP NASA 
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Figure 29H Program Output Listing : Noizexec Output (continued) 
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Figure 30 Program Output Listing : Noizexec Output (continued) 



UAAP SYSTEM FILE REQUIREMENTS 


File No. 

Logical 

Record 

Format 

Description 

Max 

Records 

5 

80 char. 

F-F-N 

Input data 

250 

6 

132 char. 

V-F-Y 

Primary output 

2500 

8 

10000 wds. 

V-U-N 

Temporary storage of 
K inverse matrix. 

” 

9 

100 wds. 

V-U-N 

Temporary storage of 
Thickness vector. 


11 

80 char. 

F-F-N 

Temporary storage of 
input data. 

250 

44 

80 char. 

F-F-N 

Temporary storage of 500 
Aero data for transfer 
to Vortex analysis 
section of code. 

50 

80 char. 

F-F-N 

Temporary storage of 
Aero and Vortex data 
for transfer to Noise 
analysis section of 
code. 

500 

51 

Ixl0**6 

wds. 

V-U-N 

Temporary file used 
in acoustic section. 


52 

Ixl0**6 

wds. 

V-U-N 

Temporart file used 
in acoustic section. 

"" 


Format : F-F-N Fixed record length with formatted records 

and no carriage control. 

V-F-Y Variable record length with formatted 

records and carriage control. 

V-u-N Variable record length with unformatted 

records and no carriage control. 

Record Length : char, characters. 

wds . words . 


Figure 31 UAAP System File Requirements 
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Figure 32 UAAP Input Data Storage Flowchart 
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PROGRAM TREE STRUCTURE 


MAIN 
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1 
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INK IN 
DATE 
CLOCK 
INPTRW 

ME WPG1 
BLDGEM 
EXIT 
BLDG 1 1 

1NPTC1 


END JOB 
EXIT 
INPTC2 
NEWPGl 
PRNTXX 
BLDCT1 

QSPLIN 

SPFUNC 

SPCOEF 

SPL5D3 

SPSVAL 

SPL5EV 

SPNVAL 

SPLSEV 

UVINT 
ME WPG] 

PRNTIN 
PR NT XX 
BDS052 
EXIT 
BDS004 

BDSOIS 

BDS014 

BDS01 0 
CAHTHA 
EXIT 
BDS01G 


EXIT 

BDS01I 

BDS020 

BDS011 

BDSOIS 

HLINE 

BDS019 

BDSOIS 


READ FROM FILE 5 AND EMCC T>C INPUT TO FILE 11 

PRINT HEADINGS /COMMANDS /DATE /TIME AT TOP OF PACE 

TOP LEVEL ROUTINE FOR BLADE GEOMETRY INPUT /CALCULATION 

INPUT DATA/CALCULATE BLADE GEOMETRY (2-D COORD OPTION) 
TREE-FIELD” LOCATION DEPENDENT INPUT ROUTINE 
TERMINATE PROGRAM — BY CALLING EXIT. 

TREE-FIELD” READER. 

PRINT ICADINGS/COHMANDS/DATE/TIME AT TOP OF PAGE 

PRINT ROUTINE FOR REAL OR INTEGER* 4 VARIABLES 

GENERATE BLADE SURFACE SHAPE A "CUT” BLADE (2-DCOORD) 

QUINT IC SPLINE INTERPOLATION 

PART OF OSPLINE INTERPOLATION 

PART OF OSPLINE INTERPOLATION 

PART OF OSPLINE IWTERPLOATION 

PART OF OSPLINE INTERPOLATION 

PART OF OSPLINE INTERPOLATION 

UNIVARIATE INTERPOLATION (NO EXTRAPOLATION) 

PRINT HEADINGS/COMMANDS/DATE /TIME AT TOP OF PAGE 
PRINT INTEGER VARIABLES 

PRINT ROUTINE FOR REAL OR INTEGER* 4 VARIABLES 
DEFINE BLADE COORD’S ON 2-D SURFACES 


DEFINE BLADE SHAPE 
DEFINE AIRFOIL SHAPE 


DETERMINE CAMBER ANGLE FROM CAMBER 


MSI 

GENERATE 

HS1 

AIRFOIL 

COORDINATES 

(DUMMY 

ROUTINE) 

EXIT 

PF1 

GENERATE 

PF1 

AIRFOIL 

COORDINATES 

(DUMMY 

ROUTINE) 

EXIT 

HS2 

GENERATE 

MS 2 

AIRFOIL 

COORDINATES 

(DUMMY 

ROUTINE) 

EXIT 

H455S 

EXIT 

PF2 

GENERATE 

PF2 

AIRFOIL 

COORDINATES 

(DUMMY 

ROUTINE) 


INTERPOLATE OVER A 4 POINT INTERVAL 
PROVIDE A PIECEWISE CUBIC SPLINE FIT OF THE DATA 
USED FOR INTERPOLATION AND EXTRAP. OF CUBIC SPLINE FIT 
CALCULATE AIRFOIL MEANLINES 

PROVIDE A PIECEWISE CUBIC SPLINE FIT OF THE DATA 
USED FOR INTERPOLATION AND EXTRAP. OF CUBIC SPLINE FIT 
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ERRSXT 

BDS025 
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BDS014 
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BDS054 

EXIT 

UNINT 

BLDGP1 

PR NT XX 

BLDGI2 

INPTC) 

END JOB 
EXIT 
INPTC2 
ME WPG! 

BLDGP1 

PRNTXX 

UNINT 

BLDG13 

INPTC1 

ENDJOB 


DETERMINE CAMBER ANCLE FROM CAMBER 

INTERPOLATE OVER A 4 POINT INTERVAL USING A VARIATION 
PART OF BLADE STACKING ROUTINES 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 


PART OF BLADE STACKING ROUTINES 
INTERPOLATE OVER A 4 POINT INTERVAL 
PART OF BLADE STACKING ROUTINES 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 


PART OF BLADE STACKING ROUTINES 

MODIFY THE AIRFOIL DATA BY CLOSING THE TRAILING/LEADING 
PRE 8 POST PROCESS INPUT I OUTPUT DATA SO THAT THE 
FIND THE LIME SEGMENTS ON UPPER S LOWER SURFACES 
FIND THE CIRCLE THAT IS TANGENT TO THREE LINES, 

MODIFY THE AIRFOIL WITH A EDGE RADIUS LARGER 
GENERATE THE POINTS THAT CLOSE THE AIRFOIL,., 

GENERATE THE POINTS THAT CLOSE THE AIRFOIL,,, 

CALCULATES ROTATIONAL COEFFICIENTS TO CONVERT 3-D 
TRANSFORMATION FROM K.Y.Z TO XS.YS.ZS 
PART OF BLADE STACKING ROUTINES 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 


PART OF BLADE STACKING ROUTINES 
PART OF BLADE STACKING ROUTINES 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 


PART OF BLADE STACKING ROUTINES 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 


UNIVARIATE INTERPOLATION OF REAL DATA 
PRINT RXYCOORD BLADE GEOMETRY INPUT 
PRINT ROUTINE FOR REAL OR INTEGER *4 VARIABLES 
INPUT BLADE GEOMETRY IN R-X-Y FORMAT (RXYCOORD) 
"TREE-FIELD” LOCATION DEPENDENT INPUT ROUTINE 
TERMINATE PROGRAM — BY CALLING EXIT. 

TREE-FIELD" READER. 

PRINT HEADINGS/COMMANDS /DATE /TIME AT TOP OF PAGE 

PRINT RXYCOORD BLADE GEOMETRY INPUT 

PRINT ROUTINE FOR REAL OR INTEGER *4 VARIABLES 

UNIVARIATE INTERPOLATION OF REAL DATA 

INPUT BLADE GEOMETRY IN X-Y-2 FORMAT (XYZCOORD) 

TREE-FIELD” LOCATION DEPENDENT INPUT ROUTINE 

TERMINATE PROGRAM — BY CALLING EXIT. 

Program Tree Structure 


Figure 33A 
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EXIT 
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ENDJOB 

EXIT 

INPTC2 

NEWPGI 


PRINT HEAD?NGS/COMKANDS/DATE/TIME AT TOP OF PAGE 
PRINT ROUTINE FOR REAL OR INTEGER" A VARIABLES 

SETUP INPUT X-V-2 COORD’S FOR BLADE CUTTING ROUTINE 

ML^LA^ EG "CUT*’ R BLAd| S FROM 3-D INPUT X-Y-2 CORRDINATES 
PART OF BLADE STACKING ROUTINES 




PART OF BLADE STACKING ROUTINES 
PART OF BLADE STACKING ROUTINES 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 


PART OF BLADE STACKING ROUTINES 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 


PRINT HEADINGS/COMMANDS/DATE/TIME AT TOP OF PAGE 
PRINT ROUTINE FOR REAL OR INTEGER" 4 VARIABLES 
UNIVARIATE INTERPOLATION OF REAL DATA 
PRINT RXYCOORD BLADE GEOMETRY INPUT 
PRINT ROUTINE FOR REAL OR INTEGER"-! VARIABLES 
DEFAULT PARAMETERS AND INPUT FOR NOISE CALC 
’TREE-FIELD" LOCATION DEPENDENT INPUT ROUTINE 
TERMINATE PROGRAM — BY CALLING EXIT. 

•TREE-FIELD" READER. _ _ 

PRINT HEADINGS /COMMANDS /DATE /TIME AT TOP OF PAGE 
PRINT ROUTINE FOR REAL OR INTEGER"* VARIABLES 

EXECUTE NOISE CALCULATION 

CARRY OUT NOISE CALCULATION - EXECUTIVE ROUTINE 
PRINT HEADINGS /COMMANDS /DATE/TIME AT TOP OF PAGE 
READS IN DATA FOR ACOUSTIC CALCULATION 
CHECK AND WRITE OUT INPUTS 

CALCULATE AUTOMATIC RADIAL INTEGRATION MESH SI2E 

MAX /MIN OF REAL ARRAY „ 

CALCULATE OMEGA INTEGRATION LIMITS FOR NEAR FIELD CALC 

I MSI SUBROUTINE 

IMSL SUBROUTINE 

IHSL SUBROUT INE 

IMSL SUBROUTINE 

MAX/MI N OF REAL ARRAY 

CALCULATE AMPLITUDE FUNCTION FOR N02KKC 

IMSL SUBROUTINE 

IMSL SUBROUTINE 

IMSL SUBROUT I 

IMSL SUBROUTINE 

CALCULATE FAR FIELD HARMONIC NOISE LEVELS 

itJTCDDni atf uoisf CALCULATION QUANTITIES TO 


PROVIDES THICKNESS SOURCE MAGNITUDE AT A GIVEN RADIUS 
CALC FOURIER TRANSFORM OF THICKNESS FOR AIRFOILS OTHER 
THICKNESS* FORMS FOR SERIES *, lb, bA, bS AND BICONVEX 
LINEAR INTERPOLATION ROUTINE 

CALCULATE FOURIER TRANSFORM OF BOUNDARY LAYER AND WAKE 
* INSTANT * QUADRUPOLE NOISE SOURCE TERM 
IMSL SUBROUTINE 

ANALYTICAL FOURIER TRANSFORM OF CHORDWISE DELTA CP 
ANAL. FOURIER TRANSFORM: REAL FUNCTION, NCP UNIFORM STEPS 
CALCULATE NEAR FIELD HARMONIC NOISE LEVELS 

LINEARLY INTERPOLATE NOISE CALCULATION QUANTITIES TO 

IMSL SUBROUTINE 
IMSL SUBROUTINE 
IMSL SUBROUTINE 
IMSL SUBROUTINE 

PROVIDES^THICKNESS SOURCE MAGNITUDE AT A GIVEN 
CALC FOURIER TRANSFORM OF THICKNESS FOR AXRFMLS OTHER 
THICKNESS FORMS FOR SERIES A, 16. 6*, bS AND BICONVEX 
LINEAR INTERPOLATION ROUTINE 

CALCULATE FOURIER TRANSFORM OF BOUNDARY LAYER AND WAKE 

•INSTANT* QUADRUPOLE NOISE SOURCE TERM 

ANALYTICAL FOURIER TRANSFORM OF CHORDWISE DELTA CP 

ANAL. FOURIER TRANSFORM: REAL FUNCTION, NCP UNIFORM STEPS 

SUM NOISE VECTORS AND PRINT FINAL RESULTS 

PRINT HEADINGS/COMMANDS/DATE/TIME AT TOP OF PAGE 

READ AND STORE DATA FOR DETERMINING FLIGHT CONDITION. 

"FREE-FIELD" LOCATION DEPENDENT INPUT ROUTINE 

TERMINATE PROGRAM — BY CALLING EXIT. 

•TREE-FIELD” READER. 

PRINT HEADINGS/COMMANDS/DATE/TIME AT TOP OF PAGE 
PRINT ROUTINE FOR REAL OR INTEGER"* VARIABLES 
TERMINATE PROGRAM — BY CALLING EXIT. 

READ AND STORE PAGE HEADERS (TITLE CARDS) 

INPUT AND DEFAULTS FOR VORTEX LIFT 
PRINT HEADINGS/COHMANDS/DATE/TIME AT TOP OF PAGE 
"FREE-FIELD" LOCATION DEPENDENT INPUT ROUTINE 
TERMINATE PROGRAM — BY CALLING EXIT. 

"FREE-FIELD" reader. 

PRINT ROUTINE FOR REAL OR INTEGER"* VARIABLES 
END OF CASE — SETUP FOR NEXT CASE. 

READ AND PRINT LSTPARMS DATA _ wli _ 

"FREE-FIELD" LOCATION DEPENDENT INPUT ROUTINE 
TERMINATE PROGRAM -- BY CALLING EXIT. 

"FREE-FIELD" READER. 

PRINT HEADINGS/COMMANDS/DATE/TIME AT TOP OF PAGE 
PRINT ROUTINE FOR REAL OR INTEGER"* VARIABLES 
PRINT COMPLEX SINGLE PRECISION VARIABLES 
INPUT AN AXSYMMETRIC FLOW FIELD. 

"FREE-FIELD" LOCATION DEPENDENT INPUT ROUTINE 
TERMINATE PROGRAM -- BY CALLING EXIT. 

"FREE-FIELD" READER. _ 

PRINT HEADINGS/COMMANDS/DATE/TIHE AT TOP OF PAGE 


Figure 33B Program Tree Structure (continued) 
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PRNTIN 

EXIT 
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DNORH 

DHNX 
1FXPLT 
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IHILU 
WAX 003 
EXIT 
HKSUB 
PI 
P2 

CHECK M 
NRK001 

JOBTIM 
NEW P<31 
DTABLP 
PRNTCD 
EXIT 
BESSAV 
EXIT 
BSJNS 
BSYS 
BLEND 

DUNLIN 

IHILU 

MUNUAR 

NRX002 

EXIT 

MRK007 

EXIT 

NRKQ04 

NRKOOS 

NRK003 

ZODEF 

DUNLIN 

NRK009 

J IKO 


FFTCXN 

FFTCB 

PLTFFO 

IHILU 


JOBTIM 

NRK006 

ZODEF 


PRINT THE AXSYWETRIC FLOW-FIELD. 

PRINT ROUTINE FOR REAL OR INTEGER "4 VARIABLES 
READ AND PRINT WAKE INPUT PARAMETERS. 

"FREE-FIELD" LOCATION DEPENDENT INPUT ROUTINE 
TERMINATE PROGRAM — BY CALLING EXIT. 

"FREE-FIELD" READER. 

PRINT HEADINGS /COMMANDS /DATE /TIME AT TOP OF PAGE 
PRINT ROUTINE FOR REAL OR INTEGER" A VARIABLES 
CONTROL PROGRAM FOR WAXE CALCULATIONS 
WAKE VELOCITY CALCULATIONS (F164WA) 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT ( DUMMY) 

PRINT HEADINCS/COMMANDS/ DATE/ TIME AT TOP OF PAGE 
PRINT WAKE INPUT PARAMETERS. 

PRINT DOUBLE PRECISION VARIABLES 
PRINT INTEGER VARIABLES 

CALCULATE F CURVE 

NORMALIZE AN ARRAY OF DOUBLE PRECISION NUMBERS 
MAX /MIN OF DOUBLE PRECISION ARRAY 
FIXED SCALE PLOTTING FOR MULTIPLE CURVES (DUMMY ROUTINE) 
IHSL SUBROUTINE 
PLOT X VS. Y (DUMMY SUBROUTINE) 

WAKE CALCULATIONS (F002WA) 

CALCULATES A CDCL VALUE AT THE GIVEN Z VALUE 
CALCULATE PART I OF THE BESSEL FUNCTON EXPRESSION 
CALCULATE PART 2 OF THE BESSEL FUNCTION EXPRESSION 
CONVERGENCE CHECKING ROUTINE FOR H SUMS 
HAKE CALCULATIONS (F1B4NR) 

RETURN THE CURRENT CPU TINE IN SEC. "1000 
PRINT HEADINGS/COMMANDS/DATE /TIME AT TOP OF PAGE 
PRINT UNIVARIATE OR BIVARIATE TABLE OF DOUBLE PREC.VAR. 
PRINT COMPLEX DOUBLE PRECISION VARIABLES 

PRECOMPUTES AND SAVES BESSEL FUNCTION ARRAYS 

IMSL SUBROUTINE 
IMSL SUBROUTINE 
BLEND AT TRAILING EDGE 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VARIABLES 
PLOT X VS. Y (DUMMY SUBROUTINE) 

COMPUTES REQUIRED VALUES OF MN AND NU 
WAKE CALCULATIONS (NEARK3) 

WAKE CALCULATIONS (NREXPS) 

CALCULATE XI BAR PART OF PSKL) (NEARA1) 

WAKE CALCULATIONS (NEARK2) 

WAKE CALCULATIONS (NEAR0) 

DEFINES THE VALUES OF 20 * DZO FOR THE BOUND TERM INTEGR. 
UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VARIABLES 
WAKE CALCULATIONS (NR INFO) 

COMPUTES BESSEL FUNCTION PRODUCT IK 
COMPUTE FAST FOURIER TRANSFORM 
IHSL SUBROUTINE 

PLOT STEADY FAST FOURIER TRANSFORMS USING IHILU SUBR . 
PLOT X VS. Y (DUMMY SUBROUTINE) 

RETURN THE CURRENT CPU TIME IN SEC. *1000 
WAKE CALCULATIONS (NEAR5 ) 

DEFINES THE VALUES OF Z 0 & DZO FOR THE BOUND TERM INTEGR. 


WAX 004 

DUNLIN 
EXIT 
NRK01 0 
EXIT 
NRK011 
NEW PCI 
1MILU 
IFXPLT 

AIRPRM 

INPTC1 

ENDJOB 

EXIT 

INPTC2 

NEWPG1 

PRNTXX 

AIRDRG 

AIRFLX 

I SOAR F 


DUNLIN 

NRK006 

1KNEW 

BSJINT 

BSJNS 

BSYINT 

BSYS 

JYNEW 

IKNWP 

JYCOMB 

FFTCXN 

FFTCB 

PLTFFO 

IHILU 

CHECK 


ISOAFL 

AIR24 

LIFT24 

L NEAR I 
DRAG24 

LNEAR1 

UNBAR 

AIROFF 

AIR24 

LIFT24 

LNEAR1 

DRAG24 

LNEAR1 

UNBAR 

CASARF 

1SOARF 

ISOAFL 

AIR24 

LIFT24 
INEAR 1 
DRAG24 
LNEAR1 
UNBAR 

AIROFF 

AIR24 
L I FT 24 
LNEAR1 


UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VARIABLES 
WAKE CALCULATIONS (NRINFM) 

COMPUTES BESSEL FUNCTION PRODUCT IK 
INTERPOLATES BESSEL FUNCTION ARRAYS 
IMSL SUBROUTINE 

INTERPOLATES BESSEL FUNCTION ARRAYS 
IMSL SUBROUTINE 

COMPUTES PRODUCT OF BESSEL FUNCTIONS JY 
COMPUTES BESSEL FUNCTION PRODUCT IK 
COMPUTES PRODUCT OF BESSEL FUNCTIONS JY 
COMPUTE FAST FOURIER TRANSFORM 
IHSL SUBROUTINE 

PLOT STEADY FAST FOURIER TRANSFORMS USING IMILU SUBR. 
PLOT X VS. Y (DUMMY SUBROUTINE) 

CONVERGENCE CHECKING ROUTINE FOR INTEGRALS 
HAKE CALCULATIONS tFOOlWA) 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VAR TABLES 
WAKE CALCULATIONS (ALNEH ) 

WAKE CALCULATIONS (HAKEBt) 

PRINT KEADIMGS/COMMANDS/DATE/TIME AT TOP OF PAGE 
PLOT X VS. Y (DUMMY SUBROUTINE) 

FIXED SCALE PLOTTING FOR MULTIPLE CURVES (DUMMY ROUTINE) 
INPUT DATA/OPT IONS FOR 2-D AIRFOIL DRAG LOOKUP 
"TREE-FIELD” LOCATION DEPENDENT INPUT ROUTINE 
TERMINATE PROGRAM -- BY CALLING EXIT. 

"FREE-FIELD” READER. 

PRINT HEADINGS/COMMANDS/DATE/TIME AT TOP OF PAGE 
PRINT ROUTINE FOR REAL OR INTEGER "4 VARIABLES 
OBTAIN 2-D DRAG FOR A GIVEN CL 

RETURN A 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B , ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 

BI -VARIATE INTERPOLATION ROUTINE 

CASCADE CORRECTION TO AIRFOIL DATA 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B >ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
B I -VAR I ATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOI L MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC. 

LINEAR INTERPOLATION AND EXTRAPOLATION 


Figure 33C Program Tree Structure (continued) 
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JOBTIM 
NEHPG1 
PRNTDP 
PR NT IN 
DTABLP 
PRNTCD 
EXIT 
DUNLIN 
BESSAV 
EXIT 
BSJNS 
BSYS 
BLEND 

DUNLIN 

IMILU 

MUNUAR 

KMATRX 

DUNLIN 

A1K2 


AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 
OBTAIN THETA-STAR FOR ALL AIRFOILS 
UNIVARIATE INTERPOLATION OF REAL DATA 
OBTAIN ANGLE OF ZERO LIFT 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR 24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI -VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B , ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR 24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

DETERMINE DELTA ALPHA ZERO LIFT FOR CASCADE CORRECT. 

B I -VAR I ATE INTERPOLATION ROUTINE 
UNIVARIATE INTERPOLATION OF REAL DATA 
CALCULATE WEINIG LIFT CURVE SLOPE FACTOR 
BI-VARIATE INTERPOLATION ROUTINE 
SWEPT AIRFOIL CORRECTION 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC. 

LINEAR INTERPOLATION AND EXTRAPOLATION 
A1R24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
Bl -VAR I ATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B , ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR 24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 

BI-VARIATE INTERPOLATION ROUTINE 

CASCADE CORRECTION TO AIRFOIL DATA 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALIE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B , ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B , ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
A1R24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI -VARIATE INTERPOLATION ROUTINE 
OBTAIN THETA-STAR FOR ALL AIRFOILS 
UNIVARIATE INTERPOLATION OF REAL DATA 
OBTAIN ANGLE OF ZERO LIFT 


APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-0 AIRFOIL LIFT FROM ALPHA , T/B , ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC. 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

DETERMINE DELTA ALPHA ZERO LIFT FOR CASCADE CORRECT. 

BI -VARIATE INTERPOLATION ROUTINE 
UNIVARIATE INTERPOLATION OF REAL DATA 
CALCULATE WEINIG LIFT CURVE SLOPE FACTOR 
BI -VARIATE INTERPOLATION ROUTINE 
LINEAR INTERPOLATION AND EXTRAPOLATION 
PRINT HEADINGS /COMMANDS/ DATE/TIME AT TOP OF PAGE 
PRINTOUT AIRFOIL OPTIONS FROM DEBUG PRINT 
PRINT INTEGER VARIABLES 

PRINT ROUTINE FOR REAL OR INTEGER" 4 VARIABLES 

PRINT ROUTINE FOR REAL OR INTEGER -4 VARIABLES 

CONTROL PROGRAM FOR AERO LOADS CALCULATIONS 
PRINT HEADINGS /COMMANDS /DATE /TIME AT TOP OF PAGE 
PRINT ROUTINE FOR REAL OR INTEGER -4 VARIABLES 
CALCULATE AERODYANMIC BLADE SWEEP. 

QUINT I C SPLINE INTERPOLATION 
PART OF QSPL1NE INTERPOLATION 
PART OF QSPLINE INTERPOLATION 
PART OF QSPLINE INTERPLOATION 

PART OF QSPLINE INTERPOLATION 

PART OF QSPLINE INTERPOLATION 
CALCULATE CGA AND FA FROM X,Y, AND Z. 

CALCULATE X,Y, AND Z FROM CGA, AND FA. 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 
CALLS ROUTINE TO GENERATE K AND WT. 

RETURN THE CURRENT CPU TIME IN SEC. "1000 

PRINT HEADINGS /COMMANDS/DATE/TIME AT TOP OF PAGE 

PRINT DOUBLE PRECISION VARIABLES 

PRINT INTEGER VARIABLES 

PRINT UNIVARIATE OR BIVARIATE TABLE OF DOUBLE PREC. VAR. 
PRINT COMPLEX DOUBLE PRECISION VARIABLES 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 
PRECOMPUTES AND SAVES BESSEL FUNCTION ARRAYS 

IMSL SUBROUTINE 
IMSL SUBROUTINE 

BLEND AT TRAILING EDGE _ 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 
PLOT X VS. Y (DUMMY SUBROUTINE) 

COMPUTES REQUIRED VALUES OF MN AND NU 

COMPUTES MATRIX OF UNSTEADY INFLUENCE C0EF ^IENTS 

inmiADiATt; I TUtTAO TUTFBPOl ATTON OF DOUBLE PREC . VAR I ABLES 


Figure 33D Program Tree Structure (continued) 
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IKSUB 

I KNEW 
BSJINT 
BSJNS 
BSYINT 
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JYNEW 
JYCOMB 
IKSUB 0 
BSYS 
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FFTCB 
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CHECK 
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MINCIN 
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BSYS 
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FFTCB 

R4ARAY 

PLTFFT 

IHILU 

CENTER 
RAD I NT 
EXIT 


RETURN THE CURRENT CPU TIME IN SEC.-1000 

SUBROUTINE TO GENERATE THE INTEGRATED MATRIX ELEMENTS 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC . VAR1 ABLES 
I MSL SUBROUTINE 
IM$L SUBROUTINE 

USED FOR IK COMBINT’N AND SLWS S * 4 IN THE Q«fl SINGLR. REG. 

CONVERGENCE CHECKING ROUTINE FOR INTEGRALS 

COMPUTE IMN INTEGRAL _ _ 

SUBROUTINE TO GENERATE THE INTEGRATED MATRIX ELEMENTS 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VAR I ABLES 

COMPUTES SUM1 IN MAKE KERNAL 

CONVERGENCE CHECKING ROUTINE FOR INTEGRALS 

COMPUTES SUM2 IN MAKE KERNAL 

COMPUTES SUMS IN MAKE KERNAL 

COMPUTE I • K* BESSEL FUNCTION USING DEBYE’S ASSYHTOTIC EXPAN. 

COMPUTE I’ * K BESSEL FUNCTION USING DEBYE’S ASSYMTOTIC EXPAN. 

COMPUTES SUM4 IN MAKE KERNAL 

USED FOR IK COMB I NT ’N AND SUMS 5 A 4 IN THE 0*0 SINGLR. REG. 
COMPUTES SUMS IN MAKE KERNAL 

USED FOR IK COMBINT’N AND SUMS 5 A 4 IN THE 0*0 SINGLR. REG. 
COMPUTES SUH6 IN MAKE KERNAL 
COMPUTES IMN INTEGRAL 

COMPUTES UNSTEADY BOUND PART OF KERNAL FUNCTION 
ROUTINE TO DEFINE THE ELEMENT IN THE MAKE KERNEL TO MHICH 
ROUTINE TO DEFINE THE ELEMENT IN TFE MAKE KERNEL TO MHICH 
UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC . VAR I ABLES 
TO CALCULATE THE W RANGE AND THE 1 1 ( MBAR , NBAR ) ARRAY 

PRECOMPUTES RADIAL INTEGRATION STATIONS 
COMPUTES M SERIES IN KERNEL 
COMPUTES IMN INTEGRAL 
COMPUTES BESSEL FUNCTION PRODUCT IK 
COMPUTES BESSEL FUNCTION PRODUCT IK 
INTERPOLATES BESSEL FUNCTION ARRAYS 
I MSL SUBROUTINE 

INTERPOLATES BESSEL FUNCTION ARRAYS 
I MSL SUBROUTINE 

COMPUTES PRODUCT OF BESSEL FUNCTIONS JY 
COMPUTES PRODUCT OF BESSEL FUNCTIONS JY 
COMPUTES BESSEL FUNCTION PRODUCT IK 
I MSL SUBROUTINE 

TO CALCLULATE THE EFFECTIVE HEIGHT OF A SLICE 
FAST FOURIER TRANSFORM 
I MSL SUBROUTINE 

TAKES 2 COMPLEX ARRAYS FROM FFT AND MAKES ♦/- REAL ARRAYS 
PLOT UNSTEADY FAST FOURIER TRANSFORMS USING IMILU SUBR . 

PLOT X VS. Y (DUMMY SUBROUTINE) 

CONVERGENCE CHECKING ROUTINE FOR INTEGRALS 

COMPUTES MING KERNEL NEAR SINGULARITY 

COMPUTES MING KERNEL NEAR SINGULARITY 

COMPUTES MING KERNEL NEAR SINGULARITY 

I MSL SUBROUTINE 

FAST FOURIER TRANSFORM 

IHSL SUBROUTINE 

TAKES 2 COMPLEX ARRAYS FROM FFT AND MAKES */- REAL ARRAYS 
PLOT UNSTEADY FAST FOURIER TRANSFORMS USING IMILU SUBR. 

PLOT X VS. Y (DUMMY SUBROUTINE) 

PERFORMS RADIAL INTEGRATION 
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CTRAP 

IMILU 

LINCG 

KMATRQ 

A1K2 
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MAKMO 
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MAKM01 

DUNLIN 

IMNSIN 

BOUND 0 

DUNLIN 

IIEXPS 

EXIT 

20DEF 

INFFNO 

I KNEW 

BSJINT 

BSJNS 

BSYINT 

BSYS 

JYNEW 

IKNMP 

JYCOMB 

IKO 

FFTCXO 

FFTCB 

PLTFFO 

IHILU 

JOBTIM 
MAKEQO 
EXIT 
DUNLIN 
BSINS 
BSKES 
CHECK 
IMNSUB 
MAKES 0 

DUNLIN 

IKPCMB 

IPKCMB 
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IKCOMB 

IMNSIN 

BOUND 

SMMAK 0 
SMMAK 
DUNLIN 
IIEXYY 
EXIT 
DEFZO 
MSUM 

INFFNS 

IKSUB 

IKNEW 

BSJINT 

BSJNS 

BSYINT 

BSYS 


LINEAR INTERPOLATION OF A UNIVARIATE CURVE 
INTEGRATE A UNIVARIATE CURVE DEFINED BY A SERIES OF 
PLOT X VS. Y (DUMMY SUBROUTINE) 

COMPUTES^MATRIX^ OF STEADY INFLUENCE COEFFICIENTS 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VARIABLES 

TO GENERATE THE INTEGRATED KATR I X ELEMENTS 

UNIVARIATE LINEAR INTERPOLATION Of DOUBLE PREC. VARIABLES 

TO GENERATE THE INTEGRATED MATRIX ELEMENTS 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 

COMPUTES IMN INTEGRAL 

COMPUTES STEADY BOUND PART OF KERNAL FUNCTION 
UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VARIABLES 
PRECOMPUTE PART OF II ASSOCIATED MITH CHM 

DEFINES THE VALUES OF ZO ( DZO FOR THE BOUND TERM INTEGR . 

COMPUTES IMN INTEGRAL 

COMPUTES BESSEL FUNCTION PRODUCT IK 

INTERPOLATES BESSEL FUNCTION ARRAYS 

I MSL SUBROUTINE 

INTERPOLATES BESSEL FUNCTION ARRAYS 
I MSL SUBROUTINE 

COMPUTES PRODUCT OF BESSEL FUNCTIONS JY 
COMPUTES BESSEL FUNCTION PRODUCT IK 
COMPUTES PRODUCT OF BESSEL FUNCTIONS JY 
COMPUTES BESSEL FUNCTION PRODUCT IK 
FAST FOURIER TRANSFORM 
I MSL SUBROUTINE 

PLOT STEADY FAST FOURIER TRANSFORMS USING IMILU SUBR. 

PLOT X VS. Y (DUMMY SUBROUTINE) 

RETURN THE CURRENT CPU TIME IN SEC. *1000 
TO GENERATE THE INTEGRATED MATRIX ELEMENTS 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 
I MSL SUBROUTINE 
I MSL SUBROUTINE 

CONVERGENCE CHECKING ROUTINE FOR INTEGRALS 
COMPUTE IMN INTEGRAL 

TO GENERATE THE INTEGRATED MATRIX ELEMENTS 
UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 
COMPUTE I • K* BESSEL FUNCTION USING DEBYE’S ASSYMTOTIC EXPAN. 
COMPUTE I * « K BESSEL FUNCTION USING DEBYE’S ASSYMTOTIC EXPAN. 
CONVERGENCE CHECKING ROUTINE FOR INTEGRALS 

USED FOR IK COMBINT’N AND SUMS 3 A 4 IN THE 0*0 SINGLR. REG. 
COMPUTES IMN INTEGRAL 

COMPUTES UNSTEADY BOUND PART OF KERNAL FUNCTION 
ROUTINE TO DEFINE THE ELEMENT IN THE WAKE KERNEL TO WHICH 
ROUTINE TO DEFINE THE ELEMENT IN THE WAKE KERNEL TO WHICH 
UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VAR 1 ABLES 


PRECOMPUTES RADIAL INTEGRATION STATIONS 
COMPUTES M SERIES IN KERNEL 
COMPUTES IMN INTEGRAL 
COMPUTES BESSEL FUNCTION PRODUCT IK 
COMPUTES BESSEL FUNCTION PRODUCT IK 
INTERPOLATES BESSEL FUNCTION ARRAYS 
I MSL SUBROUTINE 

INTERPOLATES BESSEL FUNCTION ARRAYS 
I MSL SUBROUTINE 


Figure 33E Program Tree Structure (continued) 
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FFTCB 

R4ARAY 

PLTFFT 
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CENTER 
RADI NT 
EXIT 
CUNLIN 
CTRAP 
IMILU 
L I NRG 
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PSIK2 
THICK 4 

DUNLIN 
HEX PS 
EXIT 
20DEF 
I*F 

I KNEW 

BSJINT 

BSJNS 
BSYINT 
BSYS 
JYNEW 
TIKSUB 
TIKO 
FFTCX 0 

FFTCB 

PLTFFO 

IMILU 

THICKS 

DUNLIN 
IIEX PS 
EXIT 
28DEF 
INF 

IKNEW 

BSJINT 

BSJNS 

BSYINT 

BSYS 

JYNEW 

TIKSUB 

TIKO 


COMPUTES PRODUCT OF BESSEL FACTIONS JY 
COMPUTES PRODUCT OF BESSEL FUNCTIONS JY 
COMPUTES BESSEL FUNCTION PRODUCT IK 
IMSL SUBROUTINE 

TO CALCLULATE THE EFFECTIVE HEIGHT OF A SLICE 
FAST FOURIER TRANSFORM 

IMSL SUBROUTINE _ 

TAKES 2 COMPLEX ARRAYS FROM FFT AND MAKES */- REAL ARRAYS 
PLOT UNSTEADY FAST FOURIER TRANSFORMS USING IMILU SUBR . 
PLOT X VS. Y ( DUMHY SUBROUTINE > 

CONVERGENCE CHECKING ROUTINE FOR INTEGRALS 

COMPUTES WING KERNEL NEAR SINGULARITY 

COMPUTES WING KERNEL NEAR SINGULARITY 

COMPUTES WING KERNEL NEAR SINGULARITY 

IMSL SUBROUTINE 

FAST FOURIER TRANSFORM 

IMSL SUBROUTINE 

TAKES 2 COMPLEX ARRAYS FROM FFT AND MAKES ♦/- REAL ARRAYS 
PLOT UNSTEADY FAST FOURIER TRANSFORMS USING IMILU SUBR. 
PLOT X VS. Y (DUMMY SUBROUTINE) 


PERFORMS RADIAL INTEGRATION 

LINEAR INTERPOLATION OF A UNIVARIATE CltfVE 
INTEGRATE A UNIVARIATE CURVE DEFINED BY A SERIES OF 
PLOT X VS. Y (DUMMY SUBROUTINE) 

IMSL SUBROUTINE 
COMPUTE THICKNESS VECTOR 
TO CALCULATE THE PS IV ARRAY 
COMPUTE THICKNESS VECTOR 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PR EC. VARIABLES 
PRECOMPUTE PART OF II ASSOCIATED WITH CHW 

DEFINES THE VALUES OF 20 & D20 FOR THE BOUND TERM INTEGR. 

COMPUTES IMN INTEGRAL 

COMPUTES BESSEL FUNCTION PRODUCT IK 

INTERPOLATES BESSEL FUNCTION ARRAYS 

IMSL SUBROUTINE 

INTERPOLATES BESSEL FUNCTION ARRAYS 
IMSL SUBROUTINE 

COMPUTES PRODUCT OF BESSEL FUNCT IONS JY 


FAST FOURIER TRANSFORM 
IMSL SUBROUTINE 

PLOT STEADY FAST FOURIER TRANSFORMS USING IMILU SUBR. 
PLOT X VS. V (DUMMY SUBROUTINE) 

COMPUTE THICKNESS VECTOR 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC . VARIABLES 
aocrMiDiiTP DADT ft F II ASSOCIATED WITH CHW 


DEFINES THE VALUES OF 20 * D20 FOR THE BOUND TERM INTEGR. 

COMPUTES IMN INTEGRAL 

COMPUTES BESSEL FUNCTION PRODUCT IK 

INTERPOLATES BESSEL FUNCTION ARRAYS 

IMSL SUBROUTINE 

INTERPOLATES BESSEL FUNCTION ARRAYS 
IMSL SUBROUTINE 

COMPUTES PRODUCT OF BESSEL FUNCTIONS JY 
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FFTCX 0 

FFTCB 

PLTFFO 

IMILU 

CHECK 

ANS 

LSTWC 

NEWPG1 
PRNTXX 
PRNTIN 
PRNTDP 
DUNLIN 
VEL 003 
EXIT 
NEWPG1 
PRNTDP 
PRNTXX 
PRNTIN 
BDS024 

BDS014 

ERRSXT 

BDS025 

BDS024 

BDS023 

DUNBAR 

DSLOCB 

DUNLIN 

ANSN1 

NEWPG1 

DUNLIN 

IMILU 

ANSO 

LSTWC 

NEWPG1 
PRNTXX 
PRNTIN 
PRNTDP 
DUNLIN 
VEL 003 
EXIT 
NEWPG1 
PRNTDP 
PRNTXX 
PRNTIN 
BDS024 

BDS014 

ERRSXT 

BDS025 

BDS024 

BDS023 

DUNBAR 

DSLOCB 

DUNLIN 

ANS0N1 

DUNLIN 

IMILU 

NEWPG1 

UVINT 

AJRDRG 

AIRFLX 

ISOARF 


FAST FOURIER TRANSFORM 
IMSL SUBROUTINE 

PLOT STEADY FAST FOURIER TRANSFORMS USING IMILU SUBR. 
PLOT X VS. V (DUMMY SUBROUTINE) 

CONVERGENCE CHECKING ROUTINE FOR INTEGRALS 

COMPUTE UNSTEADY POTENTIAL LIFT AND PRESSURE DISTRIB. 

INTERPOLATE THE AXSYMMETR IC FLOWFIELD FOR LST 

PRINT HEADINGS/COMMANDS/DATE/TIME AT TOP OF PAGE 

PRINT ROUTINE FOR REAL OR INTEGERS VARIABLES 

PRINT INTEGER VARIABLES 

PRINT DOUBLE PRECISION VARIABLES 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VARIABLES 
INTERPOLATE THE AXSYMMETR IC FLOWFIELD. 

PRINT HEADINGS/COMMANDS/DATE/TIME AT TOP OF PAGE 

PRINT DOUBLE PRECISION VARIABLES 

PRINT ROUTINE FOR REAL OR INTEGER 11 4 VARIABLES 

PRINT INTEGER VARIABLES 

PART OF BLADE STACKING ROUTINES 

nrockmcuT PftMPt ITAT T fMJAl FPPftR LIMIT (DUMMY) 


>ART OF BLADE STACKING ROUTINES 

)1 -VARIATE INTERPOLATION OF DOUBLE PRECISION VARIABLES 
) I FERENT 1 AT I ON OF DOUBLE PRECISION CURVE 
iniuDttTc < tkje-AD TurroDAi ATinu nr nmim F 






PRINT HEADINGS /COMMANDS/DATE/TIME AT TOP OF PAGE 
UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VARIABLES 
PLOT X VS. Y (DUMMY SUBROUTINE) 

COMPUTES STEADY POTENTIAL LIFT AND PRESSURE DISTR. 
INTERPOLATE THE AXSYMMETR I C FLOWFIELD FOR LST 
PRINT HEADINGS/COMMANDS/DATE/T I ME AT TOP OF PAGE 
PRINT ROUTINE FOR REAL OR INTEGER "4 VARIABLES 
PRINT INTEGER VARIABLES 
PRINT DOUBLE PRECISION VARIABLES 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 
INTERPOLATE THE AXSYMMETR 1C FLOWFIELD. 


PRINT HEADINGS/COMMANDS/DATE /TIME AT TOP OF PAGE 

PRINT DOUBLE PRECISION VARIABLES 

PRINT ROUTINE FOR REAL OR INTEGER "A VARIABLES 

PRINT INTEGER VARIABLES 

PART OF BLADE STACKING ROUTINES 


MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 


PART OF BLADE STACKING ROUTINES 

BI-VARIATE INTERPOLATION OF DOUBLE PRECISION VARIABLES 

D I FERENT I AT ION OF DOUBLE PRECISION CURVE 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VARIABLES 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VARIABLES 
PLOT X VS. Y (DUMMY SUBROUTINE) 

PRINT HEADINGS /COMMANDS /DATE /TIME AT TOP OF PAGE 
UNIVARIATE INTERPOLATION (NO EXTRAPOLATION) 

OBTAIN 2-D DRAG FOR A GIVEN CL 

RETURN A 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 


Figure 33F Program Tree Structure (continued) 
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AIROFF 
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UNIMT 

CLFACT 

UNBAR 

SWPARF 

ISOARF 

ISOAFL 

AIR2* 

LIFT2* 

LNEAR1 

DRAG2* 

LNEAR1 


LOOK UP CL , CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE Of LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA ,T /B»ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
B I -VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 

BI -VARIATE INTERPOLATION ROUTINE 

CASCADE CORRECTION TO AIRFOIL DATA 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL »CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR 2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC. 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
Bl -VAR I ATE INTERPOLATION ROUTINE 
OBTAIN THETA-STAR FOR ALL AIRFOILS 
UNIVARIATE INTERPOLATION OF REAL DATA 
OBTAIN ANGLE OF ZERO LIFT 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
Bl -VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR 24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
B I -VAR I ATE INTERPOLATION ROUTINE 

DETERMINE DELTA ALPHA ZERO LIFT FOR CASCADE CORRECT. 

Bl -VAR I ATE INTERPOLATION ROUTINE 
UNIVARIATE INTERPOLATION OF REAL DATA 
CALCULATE WEINIG LIFT CURVE SLOPE FACTOR 
BI-VARIATE INTERPOLATION ROUTINE 
SWEPT AIRFOIL CORRECTION 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
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NEWPG1 

AXRP02 

PRNTIN 

PRNTXX 
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CSPSUB 

BSINT 

BSJNS 

BSYS 

DUNLIN 

UVINT 

F184PE 

ERRSYT 


BI-VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC. 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 

Bl -VAR I ATE INTERPOLATION ROUTINE 

CASCADE CORRECTION TO AIRFOIL DATA 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI -VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF -AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 
OBTAIN THETA-STAR FOR ALL AIRFOILS 
UNIVARIATE INTERPOLATION OF REAL DATA 
OBTAIN ANGLE OF 2ER0 LIFT 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B , ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR 2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
B I -VAR I ATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
A1R2* DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

DETERMINE DELTA ALPHA ZERO LIFT FOR CASCADE CORRECT. 

BI-VARIATE INTERPOLATION ROUTINE 

UNIVARIATE INTERPOLATION OF REAL DATA 

CALCULATE WE IN I G LIFT CURVE SLOPE FACTOR 

BI-VARIATE INTERPOLATION ROUTINE 

LINEAR INTERPOLATION AND EXTRAPOLATION 

PRINT HEADINGS/COMMANDS/DATE/TIME AT TOP OF PAGE 

PRINTOUT AIRFOIL OPTIONS FROM DEBUG PRINT 

PRINT INTEGER VARIABLES 

PRINT ROUTINE FOR REAL OR INTEGER** VARIABLES 

PRINT ROUTINE FOR REAL OR INTEGER** VARIABLES 

COMPUTE POWER ABSORPTION DUE TO SOUND GENERATION 
INTERPOLATES BESSEL FUNCTION ARRAYS 
IMSL SUBROUTINE 
IMSL SUBROUTINE 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC . VAR I ABLES 
UNIVARIATE INTERPOLATION (NO EXTRAPOLATION) 

CALCULATE PROPELLER EFFICIENCY 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 


Figure 33G Program Tree Structure (continued) 
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ISOARF 

ISOAFL 
AIR24 

A LIFT24 

3 LNEAR1 


DRAG24 
LNEAR1 
UNBAR 
AIROFF 
AIR24 
LIFT24 
LNEAR1 
DRAG 2 4 
LNEAR1 
UNBAR 

STARC 

UN I NT 
2ER0AL 

ISOARF 

ISOAFL 

AIR24 

LIFT24 

LNEAR1 

DRAG24 

LNEARI 

UNBAR 

AIROFF 

AIR24 

LIFT24 

LNEARI 

DRAG24 

LNEARI 

UNBAR 

DZROAL 

UNBAR 

UNINT 

CLFACT 

UNBAR 

SWPARF 

ISOARF 

ISOAFL 

AIR24 

LIFT24 

LNEARI 

DRAG24 

LNEARI 

UNBAR 

AIROFF 

AIR24 

LIFT24 

LNEARI 

DRAG24 

LNEARI 

UNBAR 

CASARF 

ISOARF 

ISOAFL 

AIR24 

LIFT24 

LNEARI 

DRAG24 

LNEARI 

UNBAR 

AIROFF 

AIR24 

LIFT24 


PRINT HEADINGS /COMMANDS/ DATE/TIME AT TOP OF PAGE 
PRINT MAKE INPUT PARAMETERS. 

PRINT DOUBLE PRECISION VARIABLES 
PRINT INTEGER VARIABLES 

NORMALIZE AN ARRAY OF DOUBLE PRECISION NUMBERS 

MAX/MIN OF DOUBLE PRECISION ARRAY 

FIXED SCALE PLOTTING FOR MULTIPLE CURVES (DUMMY ROUTINE) 
IKSL SUBROUTINE 
PLOT X VS. Y (DUMMY SUBROUTINE) 

TO CALCULATE CDCLVD AT EACH OUTPUT RADIUS 

CALCULATES A CDCL VALUE AT THE GIVEN 2 VALUE 
CALCULATE PART 1 OF THE BESSEL FUNCTON EXPRESSION 
CALCULATE PART 2 OF THE BESSEL FUNCTION EXPRESSION 
CONVERGENCE CHECKING ROUTINE FOR M SIMS 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 
PERFORM SINGLE PRECISION TRAP02IDAL INTEGRATION 
PLOT X VS. Y (DUMMY SUBROUTINE) 

MOMENTUM THEORY USED TO CALC. INDUCED VEL. FOR NON-LIN. 
CORRECTS THE LINEAR MATRIX FOR MOMENTUM AXIAL INDUCED 
RETURN THE CURRENT CPU TIME IN SEC. "1000 

BLEND AT TRAILING EDGE . M 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC .VARIABLES 
PLOT X VS. Y (DUMMY SUBROUTINE) 

COMPUTES REQUIRED VALUES OF MN AND NU 
IMSL SUBROUTINE 

INCLUDES NONLINEAR EFFECTS IN L VECTOR CALCULATION,,, 

BI -VARIATE INTERPOLATION OF DOUBLE PRECISION VARIABLES 
DIFERENTI ATION OF DOUBLE PRECISION CURVE 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 
PLOT X VS. Y (DUMMY SUBROUTINE) 

PRINT HEADINGS /COMMANDS /DATE /TIME AT TOP OF PAGE 
UNIVARIATE INTERPOLATION (NO EXTRAPOLATION) 

OBTAIN 2-D DRAG FOR A GIVEN CL 

RETURN A 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC. 

LINEAR INTERPOLATION AND EXTRAPOLATION 

AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR °FF-AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC. 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 

Bl -VAR I ATE INTERPOLATION ROUTINE 

CASCADE CORRECTION TO AIRFOIL DATA 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B , ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 

AIR 24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI -VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR 0FF - AI J^ 0I L * E S^E T 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B, ETC. 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 

B I -VAR I ATE INTERPOLATION ROUTINE 

OBTAIN THETA-STAR FOR ALL AIRFOILS 

UNIVARIATE INTERPOLATION OF REAL DATA 

OBTAIN ANGLE OF ZERO LIFT „ ATA 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL»CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF -AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR 24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

DETERMINE DELTA ALPHA ZERO LIFT FOR CASCADE CORRECT. 

BI-VARIATE INTERPOLATION ROUTINE 

UNIVARIATE INTERPOLATION OF REAL DATA 

CALCULATE WEINIG LIFT CURVE SLOPE FACTOR 

BI -VAR I ATE INTERPOLATION ROUTINE 

SWEPT AIRFOIL CORRECTION „ ATA 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 

AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOX L I^SAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B , ETC. 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 

BI-VARIATE INTERPOLATION ROUTINE 

CASCADE CORRECTION TO AIRFOIL DATA 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

toOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC. 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
Bl-VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
LOOK UP * 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC. 


Figure 33H Program Tree Structure (continued) 
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LNEAR1 
DRAG 24 
LMEAR1 
UNBAR 

ST ARC 
UNI NT 
ZEROAL 
ISOARF 
1SOAFL 
AIR2A 
LIFTZft 
LMEAR1 
DRAG2A 
LMEAR1 
UNBAR 
AIROFF 
AIR24 
LIFT2A 
L NEAR I 
DRAG2A 
LNEAR 1 
UNBAR 

DZROAt 
UNBAR 
UNI NT 
CLFACT 
UNBAR 

LNEAR 1 

ME WPG 1 
AIRP02 

PR NT IN 
PRNTXX 
PRNTXX 


CSPSUB 

BSINT 

BSJNS 

BSYS 

DUNLIN 

LESWP 

VORTWT 

VTX001 

VTXINP 

EXIT 

PSPLIN 

SPFUNC 

SPCOEF 

SPL5D3 

SPSVAL 

SPL5EV 

SPHVAL 

SPL5EV 


VTXKST 

VTXDLE 

VTXCSE 

VTXCDZ 

VTXAUG 

VTXPRF 

VTXDBO 

VTXJEF 

VTXOUT 

NEHPG1 


LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR24 DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI -VARIATE INTERPOLATION ROUTINE 
OBTAIN THETA-STAR FOR ALL AIRFOILS 
UNIVARIATE INTERPOLATION OF REAL DATA 
OBTAIN ANGLE OF ZERO LIFT 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA , T/B , ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR 2ft DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
BI-VARIATE INTERPOLATION ROUTINE 

CHECK AIRFOIL DATA RETURNS FOR OFF- AIRFOIL MESSAGE 
LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC . 

LINEAR INTERPOLATION AND EXTRAPOLATION 
AIR 2ft DRAGS TABLES 

LINEAR INTERPOLATION AND EXTRAPOLATION 
B I -VAR I ATE INTERPOLATION ROUTINE 

DETERMINE DELTA ALPHA ZERO LIFT FOR CASCADE CORRECT. 

B I -VARIATE INTERPOLATION ROUTINE 
UNIVARIATE INTERPOLATION OF REAL DATA 
CALCULATE WEINIG LIFT CURVE SLOPE FACTOR 
BI -VARIATE INTERPOLATION ROUTINE 
LINEAR INTERPOLATION AND EXTRAPOLATION 
PRINT HEADINGS/COMMANDS/DATE/TIME AT TOP OF PAGE 
PRINTOUT AIRFOIL OPTIONS FROM DEBUG PRINT 
PRINT INTEGER VARIABLES 

PRINT ROUTINE FOR REAL OR INTEGER' 1 ft VARIABLES 

PRINT ROUTINE FOR REAL OR INTEGER" ft VARIABLES 

COMPUTE POWER ABSORPTION DUE TO SOUND GENERATION 
INTERPOLATES BESSEL FUNCTION ARRAYS 
1MSL SUBROUTINE 
IHSL SUBROUTINE 

UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC. VARIABLES 
CALCULATES THE LEADING EDGE SWEEP ANGLE 
WRITES OUT DATA FOR VORTEX CALCULATION,,, 

CALCULATE LOADS DUE TO VORTEX LIFT, NEW CP, CT, AND EFF. 
READ INPUT AND CALCULATE VALUES NEEDED IF AN OLD OUTPUT 

QUINTIC SPLINE INTERPOLATION 
PART OF 9SPLINE INTERPOLATION 
PART OF OSPLINE INTERPOLATION 
PART OF OSPLINE I NTERPLOAT I ON 

PART OF OSPLINE INTERPOLATION 

PART OF OSPLINE INTERPOLATION 

ADJUSTS THE K IF TOO LITTLE SWEEP, AND CALCULATES 

DISTRIBUTES THE LEADING EDGE SUCTION OVER THE CHORD 

CALCS TIP RADIAL FORCE AND SPREADS IT OVER BLADE SURFACE 

CALCULATE SPANWISE AREA WEIGHT COEFFICIENTS 

CALCULATE MAGNITUDE OF LIFT DUE TO PRESENCE OF LEADING 

RE-CALCULATE PERFORMANCE USING VORTEX LIFT CONTRIBUTION 

WRITE OUT INFORMATION FOR VORTEX CALC DEBUG 

CALCULATE EFFECTIVE ADVANCE RATIO 

WRITE OUT THE VORTEX CALCULATION RESULTS 

PRINT HEADINCS/COMMANDS/DATE/TIME AT TOP OF PAGE 


3 VTXNWT 

5 VTXWRD 

2 VTXNWT 


WRITES OUT DATA FOR ACOUSTIC CALCULATION 

WRITES OUT RADIAL LOAD DATA FOR ACOUSTIC CALCULATION 

WRITES OUT DATA FOR ACOUSTIC CALCULATION 



Figure 331 Program Tree Structure (continued) 
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CALL tS) : 

SUBROUTINE 

AIRFLX 
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CALL (S) : 
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CALL CS) : 
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CALL (S) i 
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A1RP02 
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CALL < S ) : 
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PURPOSE : 
CALL (S) : 

SUBROUTINE 

AMS 


PURPOSE t 
CALL (S) : 

SUBROUTINE 

ANSNl 


PURPOSE : 
CALL (S> r 

SUBROUTINE 

ANSO 


PURPOSE : 
CALL CS) : 

SUBROUTINE 

ANSONO 


PURPOSE : 
CALL CS) : 

SUBROUTINE 

ANSONI 


PURPOSE : 
CALL IS) : 

SUBROUTINE 

ARC TAN 


PURPOSE : 

SUBROUTINE 

A1K2 


PURPOSE : 

SUBROUTINE 

BDSA52 


PURPOSE : 
CALL <S) : 

SUBROUTINE 

BDSO 04 


PURPOSE : 
CALL (S) i 

SUBROUTINE 

BDSO 10 


PURPOSE : 
CALL (S) : 

SUBROUTINE 

BDS011 


PURPOSE : 
CALL tS) : 

SUBROUTINE 

BDS012 


PURPOSE : 

SUBROUTINE 

BDSO 1 3 


PURPOSE : 

SUBROUTINE 

BDSO 1 4 


PURPOSE 
CALL CS) : 
CALL IS) t 

SUBROUTINE 

BDS015 


PURPOSE : 
CALL CS) : 

SUBROUTINE 

BDSO 16 


PURPOSE : 
CALL tS) : 

SUBROUTINE 

BDSO 16 


PURPOSE : 

SUBROUTINE 

BDSO 19 


PURPOSE : 

SUBROUTINE 

BDS020 
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SUBROUTINE 

BDS023 
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SUBROUTINE 

BDS024 
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CALL (S) : 

SUBROUTINE 

BDS025 


PURPOSE : 
CALL (S) i 

SUBROUTINE 

6DS026 


PURPOSE : 

SUBROUTINE 

BDS052 


PURPOSE : 
CALL CS) : 

SUBROUTINE 

BDS053 


PURPOSE : 
CALL CS) : 

SUBROUTINE 

BDS0S4 


PURPOSE : 

SUBROUTINE 

BDS066 


PURPOSE : 
CALL CS) : 

SUBROUTINE 

BDS066 


PURPOSE : 






OBTAIN 2-D DRAG FOR A GIVEN CL 
AIRFLX NEWPGI AIR PC 2 PRNTXX 

RETURN A 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
ISOARF CASARF SWPARF 

CHECK AIRFOIL DATA RETURNS FOR OFF-AIRFOIL MESSAGE 
AJR24 

INPUT DAT A /OPT I ONS FOR 2-D AIRFOIL DRAG LOOKUP 
1NPTC1 NEWPG1 PRNTXX AIRDRG 
PRINTOUT AIRFOIL OPTIONS FROM DEBUG PRINT 
PRNTIN PRNTXX 

LOOK UP 2-D DRAG FOR AN INPUT VALUE OF LIFT COEFFICIENT 
LIFT24 DRAG24 

COMPUTE UNSTEADY POTENTIAL LIFT AND PRESSURE DISTRIB. 
LSTWC DUNBAR DSLOCB DUNLIN ANSNl 

NEWPGl DUNLIN IMILU 

COMPUTES STEADY POTENTIAL LIFT AND PRESSURE DISTR. 
LSTWC DUNBAR DSLOCB DUNLIN ANS0N1 NEWPG1 UVINT AIRDRG 
INCLUDES NONLINEAR EFFECTS IN L VECTOR CALCULATION, » , 
DUNBAR DSLOCB DUNLIN ANSONI NEWPG1 UVINT AIRDRG 

DUNLIN IMILU 

COMPUTE THE ARC TANGENT OF Y/X BETWEEN 0 AND 2 PI. 

CALCULATE "CUT” BLADE FROM 5-D INPUT X-Y-Z CORRDINATES 
BDS024 BDS053 BDS016 GGS252 BDS054 

BDS015 BDS104 


BDS020 

INTERPOLATE OVER A * POINT INTERVAL USING A VARIATION 
USED FOR INTERPOLATION AND EXTRAP. OF CUBIC SPLINE FIT 
DEFINE AIRFOIL SHAPE 

MLINE CAMTHA BDSO 12 BDS024 BDS020 

BDS010 BDSO 16 KS1 PF1 HS2 H455S PF2 BDS011 BDSO 19 BDSO 13 
DEFINE BLADE SHAPE 
BDSO 1 A BDS024 

: ERRSXT 

PROVIDE A PIECEWISE CUBIC SPLINE FIT OF THE DATA 

INTERPOLATE OVER A A POINT INTERVAL 

PART OF BLADE STACKING ROUTINES 

PART OF BLADE STACKING ROUTINES 

BDSO It. BDS02S 

BDSO 26 BDSQ23 

DEFINE BLADE COORD’S ON 2-D SURFACES 

EXIT BDS004 GGS107 GGS011 BDSO 2 A BDSO 5 3 BDS016 GGS252 BDS05A 

BDS02A 

MODIFY THE AIRFOIL DATA BY CLOSING THE TRAILING/LEADING 
BDS096 BDS095 BDS069 BDS086 

GENERATE THE POINTS THAT CLOSE THE AIRFOIL,., 


SUBROUTINE 

pncne? 

PURPOSE 
CALL CS) 

SUBROUTINE 

BDS095 

PURPOSE 
CALL CS) 

SUBROUTINE 

BDSO 9 6 

PURPOSE 

SUBROUTINE 

BDS0S7 

PURPOSE 

SUBROUTINE 

BDS104 

PURPOSE 
CALL CS) 

SUBROUTINE 

BESSAV 

PURPOSE 
CALL CS) 

SUBROUTINE 

BLDCT1 

PURPOSE 
CALL (S) 

SUBROUTINE 

BLDCT3 

PURPOSE 
CALL (S) 

SUBROUTINE 

BLDGEM 

PURPOSE 
CALL (S) 

SUBROUTINE 

BLDGI 1 

PURPOSE 
CALL CS) 

SUBROUTINE 

BLDGI 2 

PURPOSE 
CALL CS) 

SUBROUTINE 

BLDGI 3 

PURPOSE 
CALL CS) 

SUBROUTINE 

BLDGP1 

PURPOSE 
CALL CS) 

SUBROUTINE 

BLEND 

PURPOSE 
CALL CS) 

SUBROUTINE 

BOUND 

PURPOSE 
CALL (S) 

SUBROUTINE 

BOUNDO 

PURPOSE 
CALL (S) 

SUBROUTINE 

BSINT 

PURPOSE 
CALL (S) 

SUBROUTINE 

BSJINT 

PURPOSE 
CALL CS) 

SUBROUTINE 

BSYINT 

PURPOSE 
CALL <S> 

SUBROUTINE 

CAMTHA 

PURPOSE 
CALL IS) 

SUBROUTINE 

CASARF 

PURPOSE 
CALL (S) 

SUBROUTINE 

CENTER 

PURPOSE 

SUBROUTINE 

CHECK 

PURPOSE 

SUBROUTINE 

CHECKM 

PURPOSE 

SUBROUTINE 

CLFACT 

PURPOSE 
CALL CS) 

SUBROUTINE 

CSPSUB 

PURPOSE 
CALL (S) 

SUBROUTINE 

CTRAP 

PURPOSE 

SUBROUTINE 

CUNLIN 

PURPOSE 

SUBROUTINE 

DEFZO 

PURPOSE 

SUBROUTINE 

DMNX 

PURPOSE 

SUBROUTINE 

DNORM 

PURPOSE 
CALL CS) 

SUBROUTINE 

DRAG24 

PURPOSE 
CALL CS) 

SUBROUTINE 

DSLOCB 

PURPOSE 

SUBROUTINE 

dtablp 

PURPOSE 

SUBROUTINE 

DUNBAR 

PURPOSE 

SUBROUTINE 

DUNLIN 

PURPOSE 

SUBROUTINE 

DZROAL 

PURPOSE 
CALL (S) 


MODIFY THE AIRFOIL WITH A EDGE RADIUS LARGER 
BDS066 

FIND THE LINE SEGMENTS ON UPPER S LOWER SURFACES 
BDS097 

PRE « POST PROCESS INPUT & OUTPUT DATA SO THAT THE 
FIND THE CIRCLE THAT IS TANGENT TO THREE LINES, 

BDSQ66 

PRECOMPUTES AND SAVES BESSEL FUNCTION ARRAYS 
EXIT BSJNS BSYS 

GENERATE BLADE SURFACE SHAPE & "CUT" BLADE C2-DCOORD) 

OSPLIN UVINT NEWPG1 PRNTIN PRNTXX BDS052 EXIT UNINT 
SETUP INPUT X-Y-Z COORD’S FOR BLADE CUTTING ROUTINE 
PRNTIN BDSAS2 EXIT NEWPG1 PRNTXX UNINT 
TOP LEVEL ROUTINE FOR BLADE GEOMETRY INPUT/CALCULATION 
EXIT BLDGI1 BLDGI2 BLDGI 3 

INPUT DATA/CALCULATE BLADE GEOMETRY 12-D COORD OPTION) 

INPTC1 NEWPG1 PRNTXX BLDCT1 BLDGP1 

INPUT BLADE GEOMETRY IN R-X-V FORMAT (RXYCOORD) 

INPTC1 NEWPG1 BLDGP1 UNINT 

INPUT BLADE GEOMETRY IN X-V-Z FORMAT (XYZCOORD) 

INPTC1 NEWPGI PRNTXX EXIT BLDCT3 BLDGP1 
PRINT RXYCOORD BLADE GEOMETRY INPUT 

BLEND^AT TRAILING EDGE 
DUNLIN IMILU 

COMPUTES UNSTEADY BOUND PART OF KERNAL FUNCTION 

SMWAKO SMHAK DUNLIN 1IEXYY DEF20 MSUM WING CENTER RADINT IMILU 

COMPUTES STEADY BOUND PART OF KERNAL FUNCTION 

DUNLIN IIEXPS ZODEF 1NFFN0 

INTERPOLATES BESSEL FUNCTION ARRAYS 

BSJNS BSYS 

INTERPOLATES BESSEL FUNCTION ARRAYS 
BSJNS 

INTERPOLATES BESSEL FUNCTION ARRAYS 

BSYS 

DETERMINE CAMBER ANGLE FROM CAMBER 
EXIT 

CASCADE CORRECTION TO AIRFOIL DATA 
ISOARF ST ARC ZEROAL DZROAl CLFACT 

CONVERGENCE CHECKING ROUTINE FOR INTEGRALS 
CONVERGENCE CHECKING ROUTINE FOR M SUMS 
CALCULATE WEINIG LIFT CURVE SLOPE FACTOR 
UNBAR 

COMPUTE POWER ABSORPTION DUE TO SOUND GENERATION 
BSINT DUNLIN 

INTEGRATE A UNIVARIATE CURVE DEFINED BY A SERIES OF 

LINEAR INTERPOLATION OF A UNIVARIATE CURVE 

PRECOMPUTES RADIAL INTEGRATION STATIONS 

MAX/MIN OF DOUBLE PRECISION ARRuY 

NORMALIZE AN ARRAY OF DOUBLE PRECISION NUMBERS 

DMNX 

AIR24 DRAGS TABLES 
LNEAR1 UNBAR 

D I FERENT I AT I ON OF DOUBLE PRECISION CURVE 

PRINT UNIVARIATE OR BIVARIATE TABLE OF DOUBLE PREC.VAR. 

BI-VARIATE INTERPOLATION OF DOUBLE PRECISION VARIABLES 
UNIVARIATE LINEAR INTERPOLATION OF DOUBLE PREC . VAR I ABLES 
DETERMINE DELTA ALPHA ZERO LIFT FOR CASCADE CORRECT. 

UNBAR UNINT 


Figure 34A Subroutine References 
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SUBROUTINE 

KMATRX 


PURPOSE : 
CALL (S) : 

SUBROUTINE 

KMATR0 


PURPOSE ; 
CALL (S) : 

SUBROUTINE 

LESWP 


PURPOSE : 

SUBROUTINE 

LIFT2A 


PURPOSE : 
CALL (S) : 

SUBROUTINE 

LNEARl 


PURPOSE : 

SUBROUTINE 

LSTALS 


PURPOSE 
CALL (S> : 
CALL (S> : 

SUBROUTINE 

LSTCGA 


PURPOSE : 

SUBROUTINE 

LSTPRM 


PURPOSE l 
CALL CS) : 

SUBROUTINE 

LSTSW1 


PURPOSE : 
CALL CS) : 

SUBROUTINE 

LSTWC 


PURPOSE ; 
CALL (S> ; 

SUBROUTINE 

LSTXY2 


PURPOSE : 

SUBROUTINE 

MAIN 


PURPOSE : 
CALL CS) : 
CALL (S) : 

SUBROUTINE 

MLINE 


PURPOSE : 
CALL (S) : 

SUBROUTINE 

MNX 


PURPOSE t 

SUBROUTINE 

MOMCAL 


PURPOSE : 

SUBROUTINE 

MSUM 


PURPOSE i 
CALL (S) I 

SUBROUTINE 

MUNUAR 


PURPOSE : 

SUBROUTINE 

NEWPGI 


PURPOSE : 

SUBROUTINE 

NOISEX 


PURPOSE : 
CALL (C) : 

SUBROUTINE 

N02AIR 


PURPOSE : 

SUBROUTINE 

NQ2AMP 


PURPOSE : 
CALL (S> : 

SUBROUTINE 

NOZBLC 


PURPOSE : 

SUBROUTINE 

N02CHK 


PURPOSE 
CALL (S) 

SUBROUTINE 

N02CLC 


PURPOSE 
CALL (S) 

SUBROUTINE 

NOZDPL 


PlIR PC SE 

SUBROUTINE 

NOZFFC 


PURPOSE 
CALL (S) 

SUBROUTINE 

NOZHFL 


PURPOSE 

SUBROUT 1NE 

NOZINP 


PURPOSE 


LJ 


END OF CASE -- SETUP FOR NEXT CASE. 

TERMINATE PROGRAM -- BY CALLING EXIT. 

EXIT 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 

MACHINE DEPENDENT COMPUTATIONAL ERROR LIMIT (DUMMY) 

COMPUTE E » *X 

FAST FOURIER TRANSFORM 

FFTCB 

COMPUTE FAST FOURIER TRANSFORM 
FFTCB 

FAST FOURIER TRANSFORM 
FFTCB 

THIS ROUTINE INTERPOLATES THE F CURVE AT 2 
CALCULATE F CURVE 
DNORM IFXPLT CSINT IMILU 
TO CALCULATE CDCLVD AT EACH OUTPUT RADIUS 
EXIT HMSUB 

CALCULATE PROPELLER EFFICIENCY 

ERRSYT NEWPGI WAKPR1 EXIT FOCI F002PE DUNLIN TRAP IMILU 
CALLS ROUTINE TO GENERATE K AND NT. 

JOBTIM NEWPGI PRNTDP PRNTIN DTABLP PRNTCD EXIT DUNLIN BESSAV BLEND 

MUNUAR KMATRX KMATRO THICK ANS ANSO CSPSUB 

CORRECTS THE LINEAR MATRIX FOR MOMENTUM AXIAL INDUCED 
JOBTIM EXIT BLEND MUNUAR LINRG AMSONO CSPSUB 
TRANSFORMATION FROM X,Y,Z TO XS,YS,2S 
CALCULATES ROTATIONAL COEFFICIENTS TO CONVERT 3-D 

GGS253 GGS25A 

READ AND STORE PAGE HEADERS (TITLE CARDS) 

CALCULATES A CDCL VALUE AT THE GIVEN 2 VALUE 
PI P2 CHECKM 

GENERATE HSI AIRFOIL COORDINATES (DUMMY ROUTINE) 

EXIT 

GENERATE HS2 AIRFOIL COORDINATES (DUMMY ROUTINE) 

EXIT 

EXIT 

FIXED SCALE PLOTTING FOR MULTIPLE CURVES (DUMMY ROUTINE) 

PRECOMPUTE PART OF II ASSOCIATED WITH CHW 
EXIT 

TO CALCULATE THE W RANGE AND THE 1 1 (MBAR,NBAR ) ARRAY 
EXIT 

USED FOR IK COMB I NT * N AND SUMS 3 8 A IN THE 0*0 SINGLR. REG. 

COMPUTES BESSEL FUNCTION PRODUCT IK 
BSJINT BSYINT JYNEW 
COMPUTES BESSEL FUNCTION PRODUCT IK 
JYCOMB 

COMPUTE I ■ K ' BESSEL FUNCTION USING DEBYE’S ASSVMTOTIC EXPAN. 

COMPUTES BESSEL FUNCTION PRODUCT IK 

I KNEW JYCOMB _ 

COMPUTES BESSEL FUNCTION PRODUCT IK 
BSYS SPEC IN 

COMPUTES BESSEL FUNCTION PRODUCT IK 
PLOT X VS. Y (DUMMY SUBROUTINE) 

COMPUTES IMN INTEGRAL 

COMPUTE IMN INTEGRAL 

COMPUTES IMN INTEGRAL 

IKNEW TIKSUB TIKO FFTCXO PLTFFO 

COMPUTES IMN INTEGRAL 


IKSUB IKSUBO FFTCX R4ARAY PLTFFT 

COMPUTES IMN INTEGRAL 

IKNEW IKNWP IKO FFTCXO PLTFFO 

M FREE-F I ELD” LOCATION DEPENDENT INPUT ROUTINE 

ENDJOB INPTC2 

’’FREE -FIELD” READER. 

READ FROM FILE 5 AND EHCO THE INPUT TO FILE 11 
NEWPGI 

COMPUTE P ■ K BESSEL FUNCTION USING DEBYE'S ASSYMTOTIC EXPAN. 

LOOK UP CL, CD AS 2-D ISOLATED AIRFOIL DATA 
AIR2A AIROFF 

APPLY INPUT CORRECTION FACTORS TO AIRFOIL DATA 
ISOAFL 

RETURN THE CURRENT CPU TIME IN SEC. "1000 

COMPUTES PRODUCT OF BESSEL FUNCTIONS JY 

COMPUTES PRODUCT OF BESSEL FUNCTIONS JY 

COMPUTES MATRIX OF UNSTEADY INFLUENCE COEFFICIENTS 

DUNLIN A1K2 JOBTIM WAKE BOUND l INCG 

COMPUTES MATRIX OF STEADY INFLUENCE COEFFICIENTS 

A1K2 DUNLIN NAKMO BOUND© JOBTIM WAKEQO BOUND LINRG 

CALCULATES THE LEADING EDGE SWEEP ANGLE 

OBTAIN 2-D AIRFOIL LIFT FROM ALPHA, T/B, ETC. 

LNEAR1 

LINEAR INTERPOLATION AND EXTRAPOLATION 
CONTROL PROGRAM FOR AERO LOADS CALCULATIONS 
VORTWT VTX001 VTXNWT 

NEWPGI PRNTXX LSTSW1 DUNLIN F271 UVINT F16APE MOMCAL F271NO LESWP 
CALCULATE CGA AND FA FROM X,Y, AND 2. 

READ AND PRINT LSTPARMS DATA 
INPTC1 NEWPGI PRNTXX PRNTCS 
CALCULATE AERODYANMIC BLADE SWEEP. 

OSPLIN LSTCGA LSTXY2 

INTERPOLATE THE AXSYHMETRIC FLOWFIELD FOR 1ST 
NEWPGI PRNTXX PRNTIN PRNTDP DUNLIN VELOOS 
CALCULATE X , Y, AND 2 FROM CGA, AND FA. 

VTXPRM ENDCAS LSTPRM VELGRD WAKPRM WAX 001 AIRPRM LSTALS 

INK IN DATE CLOCK INPTRW BLDGEM N02PRM NOISEX RUNPRM ENDJOB HEADER 

CALCULATE AIRFOIL MEANLINES 

BDS019 BDSO 1 3 

MAX/MIN OF REAL ARRAY 

MOMENTUM THEORY USED TO CALC. INDUCED VEL . FOR NON- LIN. 

COMPUTES M SERIES IN KERNEL 
IWFFNS CHECK 

COMPUTES REQUIRED VALUES OF MN AND NU 

PRINT HEADINGS/COMMANDS/DATE/TIME AT TOP OF PAGE 

EXECUTE NOISE CALCULATION 

N02CLC 

FIND AIRFOIL TYPE AT 2 

CALCULATE AMPLITUDE FUNCTION FOR NOZKKC 
BS1NS BSKES WSJNS BSYS 

CALCULATE FOURIER TRANSFORM OF BOUNDARY LAYER AND WAKE 
CHECK AND WRITE OUT INPUTS 
N02NJ2 NOZKKC 

CARRY OUT NOISE CALCULATION - EXECUTIVE ROUTINE 
NEWPGI N02INP N02CHK N02FFC N02NFC N02SMS 
ANALYTICAL FOURIER TRANSFORM OF CMORDW1SE DELTA CP 
CALCULATE FAR FIELD HARMONIC NOISE LEVELS 
NOZTRP NOZAIR NOZHPL NOZBLC NOZQPL BSJNS N02DPL NOZRDL 
THICKNESS FORMS FOR SERIES A, lb, b6, b5 AND BICONVEX 
READS IN DATA FOR ACOUSTIC CALCULATION 


Figure 34B Subroutine References (continued) 



SUBROUTINE N02KKC 
SUBROUTINE N02HPL 


SUBROUTINE N02NFC 


SUBROUTINE 

SUBROUTINE 

SUBROUTINE 


SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 


SUBROUTINE 


SUBROUTINE 

SUBROUTINE 


SUBROUTINE 

SUBROUTINE 


I SUBROUTINE I 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 

SUBROUTINE 

SUBROUTINE 

SUBROUTINE 


PURPOSE : CALCULATE OMEGA INTEGRATION LIMITS FOR NEAR FIELD CALC 
CALL (S) : BSINS BSKES BSJNS BSVS MNX N02AMP 
PURPOSE : PROVIDES THICKNESS SOURCE MAGNITUDE AT A GIVEN RADIUS 
CALL tS) : N02NHF 

PURPOSE : CALCULATE NEAR FIELD HARMONIC NOISE LEVELS 

CALL <S) : EXI? F>L N02TRP BSINS BSKES BSJNS BSYS N02AIR N02MPL N02BLC N020PL 
PURPOSE : CALC FOURIER TRANSFORM OF THICKNESS FOR AIRFOILS OTHER 
CALL (S) : N02HFL UNlIN 

PURPOSE : CALCULATE AUTOMATIC RADIAL INTEGRATION MESH SUE 
CALL (S) : MNX 

PURPOSE : DEFAULT PARAMETERS AND INPUT FOR NOISE CALC 
CALL CS) : INPTC1 NEWPG1 PRNTXX 

PURPOSE : 'INSTANT' QUADRUPOLE NOISE SOURCE TERM 

PURPOSE : ANAL. FOURIER TRANSFORM: REAL FUNCTION, NCP UNIFORM STEPS 
PURPOSE : SUM NOISE VECTORS AND PRINT FINAL RESULTS 
CALL (S) : NEWPG1 

PURPOSE : LINEARLY INTERPOLATE NOISE CALCULATION QUANTITIES TO 
CALL (S) : EXIT 

PURPOSE : MAKE CALCULATIONS CF1BANR) 

CALL (S) : JOBTIH NEWPG1 DTABLP PRNTCD EXIT BESSAV BLEND MUNUAR NRK002 
PURPOSE : WAKE CALCULATIONS CNEARK3) 

CALL IS) : EXIT NRK007 NRKOOA NRK005 MRK003 JOBTIM NRKOOfc 
PURPOSE : WAKE CALCULATIONS (NEARO) 

CALL (S) : 20DEF DUNLIN NRK009 

PURPOSE : CALCULATE XI BAR PART OF PSICL) (NEARA1) 

PURPOSE : WAKE CALCULATIONS ( WEARK2) 

PURPOSE : WAKE CALCULATIONS < NEAR3 ) 

CALL (S) : 20DEF DUNLIN NRK008 CHECK 
PURPOSE : WAKE CALCULATIONS CNREXPS) 

CALL (S) : EXIT 

PURPOSE t WAKE CALCULATIONS (NRINFM) 

CALL (S) : IKNEW IKNWP FFTCXN PLTFFO 
PURPOSE : HAKE CALCULATIONS (NRIWFO) 

CALL (S) : IK 0 FFTCXN PLTFFC 
PURPOSE : WAKE CALCULATIONS (ALNEW ) 

CALL (S) : EXIT NRK011 NEWPG1 IMILU IFXPLT 
PURPOSE : WAKE CALCULATIONS IWAKE66) 

PURPOSE : GENERATE PF1 AIRFOIL COORDINATES (DUMMY ROUTINE) 

PURPOSE 5 : GENERATE PF2 AIRFOIL COORDINATES (DUMMY ROUTINE) 

PURPOSE 1 : PLOT UNSTEADY FAST FOURIER TRANSFORMS USING IMILU SUBR . 

PURPOSE 1 * PLOT U STEADY FAST FOURIER TRANSFORMS USING IMILU SUBR. 

CALL (S> : IMILU 

PURPOSE : PRINT COMPLEX DOUBLE PRECISION VARIABLES 

PURPOSE : PRINT COMPLEX SINGLE PRECISION VARIABLES 

PURPOSE : PRINT DOUBLE PRECISION VARIABLES 
PURPOSE : PRINT INTEGER VARIABLES 

PURPOSE : PRINT ROUTINE FOR REAL OR INTEGER* A VARIABLES 
PURPOSE : TO CALCULATE THE PSIV ARRAY 

PURPOSE : CALCULATE PART 1 OF THE BESSEL FUNCTON EXPRESSION 

PURPOSE : CALCULATE PART 2 OF THE BESSEL FUNCTION EXPRESSION 

PURPOSE : QUINT I C SPLINE INTERPOLATION 

CALL (S> : SPFUNC SPCOEF SPSVAL SPMVAL 

PURPOSE : PERFORMS RADIAL INTEGRATION 

CALL (S) : EXIT CUNL IN CTRAP 

PURPOSE : FUNCTION TO COMPUTE RJ AND DERIVATES 

PURPOSE : READ AND STORE DATA FOR DETERMINING FLIGHT CONDITION. 

CALL CS) : INPTC1 NEWPG1 PRNTXX 


SUBROUTINE 

R6ARAY 

PURPOSE 

SUBROUTINE 

SMWAK 

PURPOSE 

SUBROUTINE 

SMWAKO 

PURPOSE 

SUBROUTINE 

SPCOEF 

PURPOSE 
CALL CS) 

SUBROUTINE 

SPECIN 

PURPOSE 

SUBROUTINE 

SPFUNC * 

PURPOSE 

SUBROUTINE 

SPL5D3 

PURPOSE 

SUBROUTINE 

SPL5EV 

PURPOSE 

SUBROUTINE 

SPMVAL 

PURPOSE 
CALL CS) 

SUBROUTINE 

SPSVAL 

PURPOSE 
CALL CS) 

SUBROUTINE 

STARC 

PURPOSE 
CALL <S) 

SUBROUTINE 

SUM1 

PURPOSE 

SUBROUTINE 

SUM2 

PURPOSE 

SUBROUTINE 

SUM 3 

PURPOSE 
CALL CS) 

SUBROUTINE 

SUM 

PURPOSE 
CALL CS) 

SUBROUTINE 

SUMS 

PURPOSE 
CALL CS) 

SUBROUTINE 

SUM 6 

PURPOSE 

SUBROUTINE 

SWPARF 

PURPOSE 
CALL CS) 

SUBROUTINE 

THICK 

PURPOSE 
fTALL CS) 

SUBROUTINE 

THICKO 

PURPOSE 
CALL CS) 

SUBROUTINE 

THICKS 

PURPOSE 
CALL CS) 

SUBROUTINE 

TIKSUB 

PURPOSE 

SUBROUTINE 

TIKO 

PURPOSE 

SUBROUTINE 

TRAP 

PURPOSE 

SUBROUTINE 

UNBAR 

PURPOSE 

SUBROUTINE 

UNI NT 

PURPOSE 

SUBROUTINE 

UNLIN 

PURPOSE 

SUBROUTINE 

UVINT 

PURPOSE 

SUBROUTINE 

VELGRD 

PURPOSE 
CALL tS) 

SUBROUTINE 

VEL 002 

PURPOSE 
CALL CS) 

SUBROUTINE 

VEL003 

PURPOSE 
CAUL CS) 

SUBROUTINE 

VORTEX 

PURPOSE 
CALL (S) 

SUBROUTINE 

VORTWT 

PURPOSE 

SUBROUTINE 

VTXAUG 

PURPOSE 

SUBROUTINE 

VTXCD2 

PURPOSE 

SUBROUTINE 

VTXCSE 

PURPOSE 

SUBROUTINE 

VTXDBO 

PURPOSE 

SUBROUTINE 

VTXDLE 

PURPOSE 

SUBROUTINE 

VTXINP 

PURPOSE 
CALL (S) 

SUBROUTINE 

VI XJEF 

PURPOSE 

SUBROUTINE 

VTXKST 

PUR POSE 

SUBROUTINE 

VTXNWT 

PURPOSE 

SUBROUTINE 

VTXOUT 

PURPOSE 
CALL CS) 

SUBROUTINE 

VTXPRF 

PURPOSE 


TAKES 2 COMPLEX ARRAYS FROM FFT AND MAKES ♦/- REAL ARRAYS 
ROUTINE TO DEFINE THE ELEMENT IN THE WAKE KERNEL TO WHICH 
ROUTINE TO DEFINE THE ELEMENT IN THE WAKE KERNEL TO WHICH 
PART OF QSPLINE INTERPOLATION 
SPL5D3 

TO CALCLULATE THE EFFECTIVE HEIGHT OF A SLICE 
PART OF QSPLINE INTERPOLATION 
PART OF QSPLINE INTER PL OAT I ON 
PART OF QSPLINE INTERPOLATION 


CALL (S) : SPL5EV 

PURPOSE : OBTAIN THETA-STAR FOR ALL AIRFOILS 
CALL <S) J UNI NT 

PURPOSE : COMPUTES SUM1 IN WAKE KERNAL 

PURPOSE : COMPUTES SUM2 IN WAKE KERNAL 

PURPOSE : COMPUTES SUM3 IN WAKE KERNAL * 

CALL (S) : IKPCMB IPKCMB 

PURPOSE : COMPUTES SUM4 IN WAKE KERNAL 

CALL <S) : IK COMB 

PURPOSE : COMPUTES SUMS IN WAKE KERNAL 
CALL CS) : IKCOMB 

PURPOSE : COMPUTES SUMfc IN WAKE KERNAL 

PURPOSE : SWEPT AIRFOIL CORRECTION 

CALL CS) : ISOARF CASARF LNEAR1 

PURPOSE : COMPUTE THICKNESS VECTOR 

fT ALL CS) : PSIK2 THICKO THICKS 

PURPOSE : COMPUTE THICKNESS VECTOR 

CALL (S) : DUNLIN IIEXPS 20DEF INF 

PURPOSE : COMPUTE THICKNESS VECTOR 

CALL CS) t DUNLIN IIEXPS 20DEF INF CHECK 

PURPOSE : 

PURPOSE : PERFORM SINGLE PRECISION TRAP02IDAL INTEGRATION 

PURPOSE : B1 -VARIATE INTERPOLATION ROUTINE 

PURPOSE : UNIVARIATE INTERPOLATION OF REAL DATA 

PURPOSE : LINEAR INTERPOLATION ROUTINE 

PURPOSE : UNIVARIATE INTERPOLATION (NO EXTRAPOLATION) 

PURPOSE ; INPUT AN AXSYMMETR I C FLOW FIELD. 

CALL (S) : INPTC1 NEWPG1 VEL002 

PURPOSE : PRINT THE AXSYMMETR IC FLOW-FIELD. 

CALL CS) : PRNTXX 

PURPOSE : INTERPOLATE THE AXSYMMETR I C FLOWF1ELD. 

CAUL (S) : EXIT NEWPG1 PRNTDP PRNTXX PR NT IN BDS02A 
PURPOSE : EXECUTE VORTEX CALCULATION 
CALL (S) : VTX001 

PURPOSE : WRITES OUT DATA FOR VORTEX CALCULATION,,, 

PURPOSE : CALCULATE MAGNITUDE OF LIFT DUE TO PRESENCE OF LEADING 

PURPOSE : CALCULATE SPANWISE AREA WEIGHT COEFFICIENTS 

PURPOSE : CALCS TIP RADIAL FORCE AND SPREADS IT OVER BLADE SURFACE 

PURPOSE : WRITE OUT INFORMATION FOR VORTEX CALC DEBUG 

PURPOSE : DISTRIBUTES THE LEADING EDGE SUCTION OVER THE CHORD 

PURPOSE : READ INPUT AND CALCULATE VALUES NEEDED IF AN OLD OUTPUT 

CALL (S) : EXIT QSPLIN 

PURPOSE ! CALCULATE EFFECTIVE ADVANCE RATIO 

PURPOSE : ADJUSTS THE K IF TOO LITTLE SWEEP, AND CALCULATES 
PURPOSE : WRITES OUT DATA FOR ACOUSTIC CALCULATION 
PURPOSE ! WRITE OUT THE VORTEX CALCULATION RESULTS 

PURPOSE 1 : RE-CALCULATE PERFORMANCE USING VORTEX LIFT CONTRIBUTION 

Figure 34C Subroutine References (continued) 
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SUBROUTINE 

VT X PR M 

PURPOSE 
CALL i S) 

SUBROUTINE 

VTXWRD 

PURPOSE 

SUBROUTINE 

V7X001 

PURPOSE 
CALL (S) 
CALI (S) 

SUBROUTINE 

WAKE 

PURPOCE 
CALL tS> 

SUBROUTINE 

WAKEQO 

PUR POSE 
CALL CS> 

SUBROUTINE 

WAKESN 

PURPOSE 
CALL (S) 

SUBROUTINE 

WAKES 0 

PURPOSE 
CALL (S) 

SUBROUT I WE 

WAK MO 

PURPOSE 
CALL tS) 

SUBROUTINE 

WAKMG1 

PURPOSE 
CALL (S) 

SUBROUTINE 

WAKPRM 

PURPOSE 
CALL <S) 

SUBROUTINE 

WAK PR 1 

PURPOSE 
CALL <S) 

SUBROUTINE 

WAK001 

PURPOSE 
CALL tS) 

SUBROUTINE 

WAX 002 

PURPOSE 
CALL (S) 

SUBROUTINE 

WAK 003 

PURPOSE 
CALL (S) 

SUBROUTINE 

WAK 004 

PURPOSE 
CALL tS) 

SUBROUTINE 

WING 

PURPOSE 
CALL CS) 

SUBROUTINE 

WINGF 

PURPOSE 
CALL <S) 

SUBROUTINE 

WINGIN 

PURPOSE 
CALL <S) 

SUBROUTINE 

ZEROAL 

PURPOSE 
CALL tS) 

SUBROUTINE 

Z0DEF 

PURPOSE 


INPUT AND DEFAULTS FOR VORTEX LIFT 
NEWPGl 1NPTC1 PRNTXX 

WRITES OUT RADIAL LOAD DATA FOR ACOUSTIC CALCULATION 
CALCULATE LOADS DUE TO VORTEX LIFT, NEW CP, CT , AND EFF . 

VTXNWT VTXWRD 

VTXINP VTXKST VTXDLE VTXCSE VTXCDZ VTXAUG VTXPRF VTXDBO VTKJEF VTXOUT 

SUBROUTINE TO GENERATE THE INTEGRATED MATRIX ELEMENTS 

EXIT DUNLIN BSINS BSKES IK COMB CHECK IMNSUB WAKESN 

TO GENERATE THE INTEGRATED MATRIX ELEMENTS 

EXIT DUNLIN BSINS BSKES CHECK IMNSUB WAKESO 

SUBROUTINE TO GENERATE THE INTEGRATED MATRIX ELEMENTS 

DUNLIN SUM1 CHECK SUM2 SUMJ SUMA SUMS SUM6 IMNSIN 

TO GENERATE THE INTEGRATED MATRIX ELEMENTS 

DUNLIN IKPCMB IPKCMB CHECK IKCOMB IMNSIN 

TO GENERATE THE INTEGRATED MATRIX ELEMENTS 

DUNLIN WAKM01 

TO GENERATE THE INTEGRATED MATRIX ELEMENTS 
DUNLIN IMNSIN 

READ AND PRINT WAKE INPUT PARAMETERS. 

1NPTCI NEWPG1 PRNTXX 
PRINT WAKE INPUT PARAMETERS. 

PRNTDP PRNTIN 

CONTROL PROGRAM FOR WAKE CALCULATIONS 
WAK 0 0 2 

WAKE VELOCITY CALCULATIONS (F1B4WA) 

ERRSYT NEWPG1 WAKPR1 EXIT F001 WAK00S 
WAKE CALCULATIONS (F002WA) 

EXIT HMSUB NRK001 WAK004 NRK010 

WAKE CALCULATIONS (F001WA) 

DUNLIN EXIT 

COMPUTES WING KERNEL NEAR SINGULARITY 
WINGIN 

COMPUTES WING KERNEL NEAR SINGULARITY 
BSYS 

COMPUTES WING KERNEL NEAR SINGULARITY 

WINGF FFTCX R4ARAY PlTFFT 

OBTAIN ANGLE OF ZERO LIFT 

ISOARF 

DEFINES THE VALUES OF 20 & DZO FOR THE BOUND TERM INTEGR. 


! 
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Figure 34D 
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SUBROUTINE CROF* REFERENCE LISTING 
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SUM. SUBROUTINE NAMES 


AIRDRG IS CALLED BY : AIRPRH AMSO ANSONO 

AIRFLX IS CALLED BY : AIRDRG 

AIROFF IS CALLED BY : ISOAFL 

AIRPRH IS CALLED BY : MAIN 

AIRP02 IS CALLED BY * AIRDRG 

AIR24 IS CALLED BY : AIROFF ISOAFL 

AMS IS CALLED BY : F271 

AMSNl IS CALLED BY t ANS 

AMSO IS CALLED BY : F27I 

AMS DMO IS CALLED BY : F271NO 

AMS ONI IS CALLED BY t AMSO ANSONO 

ARCTAN IS CALLED BY : 

AIK 2 IS CALLED BY * KNATRX KHATRO 

BDSA52 IS CALLED BY r BLDCTS 

BDS004 IS CALLED BY t BDSDS2 

BDSOI 0 IS CALLED BY t BDSOU 

BDS011 IS CALLED BY t BDSOU 

BDSOI 2 IS CALLED BY t BDSOU 

BDSOI S IS CALLED BY t BDSOU KLINE 

BDSOU IS CALLED BY ! BDSOU 

BDSOI 5 IS CALLED BY l BDS004 

BDSOU IS CALLED BY t BDSAS2 BDS052 BDS024 

BDSOU IS CALLED BY i BDSOU 

BDSOU IS CALLED BY i BDSOU MLINE 

BDS020 IS CALLED BY i BDSOU BDSOU 

BDS02S IS CALLED BY : BDS025 

BDS024 IS CALLED BY I VEL005 BDSA52 BDSOU BDSOIS BDS052 BDSOS 3 

BDS025 IS CALLED BY : BDS024 

BDS020 IS CALLED BY t BDS025 

BDS052 IS CALLED BY i BLDCTl 

BDS0S3 IS CALLED BY t BDSAS2 BDS052 

BDSQ54 IS CALLED BY t BDSA52 BDS052 

BDSOBt IS CALLED BY t BDS104 

BOSOM IS CALLED BY t BDSOBt BDS060 

BDS0G9 IS CALLED BY t BDSOBt 

BDS095 IS CALLED BY t BDSOBt 

BDS09t IS CALLED BY t BDSOAt 

BDS097 IS CALLED BY l BDS095 

BDS104 IS CALLED BY I BDS004 

BESSAV IS CALLED BY : F271 NRK001 

BLDCTl IS CALLED BY s BLDGI1 

BLDCT3 IS CALLED BY : BLDG IS 

BLDGEH IS CALLED BY : MAIN 

BLDG 11 IS CALLED BY : BLDGEH 

BLDGI2 IS CALLED BY * BLDGEH 

BLDGIS IS CALLED BY j BLDGEH 

BLDGP1 IS CALLED BY J BLDGI1 BLDGI2 BLDGIS 

BLEND IS CALLED BY t F271 F271NO NRK001 

BOUND IS CALLED BY t KHATRX KHATRO 

BOUND 0 IS CALLED BY : KHATRO 

BSIMT IS CALLED BY : CSPSUB 

BSJ1NT IS CALLED BY t I KNEW 

BSYIWT IS CALLED BY * 1KNEH 

CAMTHA IS CALLED BY : BDSOU BDSOU 

CASARF IS CALLED BY J AIRFLX SHPARF 

CENTER IS CALLED BY : BOUND 


IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
, IS CALLED BY 
i IS CALLED BY 
> IS CALLED BY 
I IS CALLED BY 
| IS CALLED BY 


IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
I IS CALLED BY 
i IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
; IS CALLED BY 
: is CALLED BY 
IS CALLED BY 
> IS CALLED BY 
, IS CALLED BY 
r IS CALLED BY 
l IS CALLED BY 
l IS CALLED BY 
i IS CALLED BY 
l IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
r IS CALLED BY 
i IS CALLED BY 
i IS CALLED BY 
* IS CALLED by 

is called by 

IS CALLED BY 
B IS CALLED BY 
IS CALLED BY 
0 IS CALLED BY 
IS CALLED BY 
IS CALLED BY 
N IS CALLED BY 
B IS CALLED BY 
IS CALLED BY 
S IS CALLED BY 

0 IS CALLED BY 

1 IS CALLED BY 

2 IS CALLED BY 
H IS CALLED BY 
IB IS CALLED BY 


THICXS HAKEQO HAKES 0 HAKE HAKESN MSUM 


NRKOOt 

msus 

CASARF 

F271 

RAD I NT 

RAD I NT 

BOUND 

DMORH 

F001 

AIR24 

AMSO 

F271 

AMSO 

HAKES 0 

KHATRX 

lstwc 

CASARF 

MAIN 

MAIN 

BDSOU 

F1B4PE 
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MAXHO IS CALLED 
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Subroutine Cross References (continued) 
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